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ABSTRACT 


The oblique wing concept for transonic aircraft has been 
proposed to reduce drag. This work investigates the dynamic stability 
of the aircraft by analytically determining the stability derivatives 
at angles of skew ranging from 0° and 45° and using these stability 
derivatives in a linear analysis of the coupled aircraft behavior. The 
stability derivatives were obtained using a lifting line aerodynamic 
theory and found to give reasonable agreement with derivatives developed 
in a Boeing study for the same aircraft. 

In the dynamic analysis, no instability or large changes 

o o 

occurred in the root locations for skew angles varying from 0 to 45 
with the exception of roll convergence. The damping in roll, however, 
decreased by an order of magnitude. Furthermore, rolling was a promi- 
nent feature of all the oscillatory .mode shapes at high skew angles. 


i 



TABLE OF CONTENTS 


Page 

ABSTRACT . i 

LIST OF SYMBOLS vi 

LIST OF FIGURES xi 

LIST OF TABLES xvi 

I. INTRODUCTION 1 

1.1 Background ...... \ 

1.2 Synopsis and Contributions 3 

II. AERODYNAMICS • • 5 

2.1 Generalities 5 

2.2 Reference Axes • 8 

2.2.1 Body Axes (x^ , y^ , z^) 8 

2.2.2 Stability Axes (x g , y g , z^) 9 

2.3 Strip Theory 9 

2.3.1 Variations in Aerodynamic Forces due to a 

Perturbation OL 11 

2. 3. 1.1 Section Lift Slope • • 11 

2.3. 1.2 Section Downwash 12 

2.3. 1.3 Section Induced Drag . 14 

2.3.2 Evaluation of Stability Derivatives by 

Means of Strip Theory ..... 18 

2.4 Lifting Line Theory 19 

2.4.1 Introduction 19 

2.4.2 Steady State Loading on an Airplane with 

an Oblique Wing 23 

ii 



Page 

2.4.3 Section Final Angle of Attack . 26 

2.4.4 Section Induced Drag . 28 

2.5 An Empirical Correction to Schrenk's Method ... 29 

2.6 Summary - 35 

III. STABILITY DERIVATIVES 36 

3.1 Basic Features 36 

3.1.2 Side Force 39 

3.1.3 Wing Rotation 49 

3.2 Derivatives 54 

3.2.1 General Methodology 54 

3.2.2 Oi , p , q Derivatives ........... 56 

3.2.3 [5 Derivatives 61 

3.2.4 r Derivatives 62 

3.3 Summary 65 

IV. EQUATIONS OF MOTION 66 

4.1 Rigid Body 66 

4.2 Effects of Flexibility 71 

V. NUMERICAL ANALYSIS 75 

5.1 Model Description 75 

5.2 Aerodynamic Results 78 

5.2.1 Comparison Between Strip Theory and Linear 

Theory (Rigid Case) ..... 78 

5.2.2 Comparison Between Rigid and Elastic Wing 

(Linear Theory) 88 

5.3 Stability Derivatives and their Influence on the 

Natural Modes 120 

iii 



Page 

5.4 Stability Derivatives ........ 126 

5.5 Dynamic Results . 150 

5.5.1 Natural Modes . . 150 

5.5.2 Transient Response * . 161 

5.6 Summary 164 

VI. CONCLUSIONS 165 

APPENDIX A . . . 166 

A. l Relations Between Dimensional and Nondimensional . . 

Derivatives (Stability Axes) ... 166 

APPENDIX B 168 

B. l u Derivatives 168 

B.2 {J Derivatives 168 

B.3 CH Derivatives 174 

B.4 p Derivatives 179 

B.5 q Derivatives 183 

B. 6 r Derivatives . 186 

APPENDIX C 189 

C. l Similarity Transformation Between Inertia Matrices 189 

APPENDIX D 191 

D. l Downwash Matrix [S^3 191 

D.2 Evaluation of K„ Elements 196 

D.3 Structures Fundamentals 205 

D.4 Wing Twist 223 

D.4.1 Twist due to Control Surface Deflection . . 224 

iv 



Page 


D. 4.2 Apparent Twist due to Rolling and Pitching 

Velocities . 229 

D. 5 Un symmetrical Flight Conditions . 231 

REFERENCES 240 


v 



LIST OF SYMBOLS 


0 

m. 

b' 

b o 

c or C 


'R, 


a 


m 


m 


m 


m 


m 


a 


m 


'a 


2-D Section Lift Slope 

Aspect Ratio 

Actual Wing Span 

Unyawed Wing Span 

Local Wing Chord 

Pivot - C.G. distance (Fig. 3.5) 

Mean Aerodynamic Chord for Unyawed Wing 
Root Chord for Unyawed Wing 

Rolling Moment Stability Derivative due to Roll Rate 
Rolling Moment Stability Derivative due to Pitch Rate 
Rolling Moment Stability Derivative due to Yaw Rate 


Rolling Moment Stability Derivative due to Angle of Attack 
Rolling Moment Stability Derivative due to Sideslip 
Pitching Moment Coefficient 


Pitching Moment Stability Derivative due to Roll Rate 

Pitching Moment Stability Derivative due to Pitch Rate 

Pitching Moment Stability Derivative due to Yaw Rate 

Pitching Moment Stability Derivative due to Angle of Attack 

Pitching Moment Stability Derivative due to Angle of Attack Rate 

Pitching Moment Stability Derivative due to Sideslip 


vi 




Yawing Moment Stability Derivative due to Roll Rate 
Yawing Moment Stability Derivative due to Pitch Rate 
Yawing Moment Stability Derivative due to Yaw Rate 
Yawing Moment Stability Derivative due to Angle of Attack 
Yawing Moment Stability Derivative due to Sideslip 
Horizontal Force Coefficient 

Side Force Stability Derivative due to Angle of Attack 
Side Force Stability Derivative due to Roll Rate 
Side Force Stability Derivative due to yaw rate 
Side Force Stability Derivative due to Sideslip 
Vertical Force Stability Derivative due to Roll Rate 
Vertical Force Stability Derivative due to Pitch Rate 


q 


f C L ( ^ ] B 

[c L (y)] F 

C T (y) 

[C. L (y)]Tot 


Vertical Force Stability Derivative due to 

Vertical Force Stability Derivative due to 

Vertical Force Stability Derivative due to 

Lift Coefficient 
Section Basic-Lift Coefficient 
Flat Wing Section Lift Coefficient 
Section Lift Slope 

Section Total Lift Coefficient [Ref. Eq. 
Induced Drag Coefficient 


Angle of Attack 
Angle of Attack Rate 
Sideslip 


2 . 5 ] 


1 


vii 



C D (y) Section Induced Drag Coefficient 

(C_ (y) ] Section Induced Drag due to Basic Lift 
D f B 


D 

D. 


<Vy 


d. 

1 

FWLD 

F 


(F ) 
v y w 

(P A 

I 

I 

XX 

I' 

XX 


yy 

i’ 

yy 


zz 

l ! 

zz 


xy 

I' 

xy 


xz 

I* 

xz 


Drag Force 

Section Induced Drag Force 

Side Component of the Section Induced Drag Force 

Section Drag 

Section Induced Drag 

Flat Wing Lift Distribution 

External Forces 

x-Component of the Aerodynamic Force 
y-Component of the Aerodynamic Force 
z-Component of the Aerodynamic Force 
Side Force due to Induced Drag 
Side Force due to Wing Lift 
Inertia Matrix 

Roll Moment of Inertia in Body Axes 
Roll Moment of Inertia in Stability Axes 
Pitch Moment of Inertia in Body Axes 
Pitch Moment of Inertia in Stability Axes 
Yaw Moment of Inertia in Body Axes 
Yaw Moment of Inertia in Stability Axes 
Roll-Pitch Product of Inertia in Body Axes 
Roll-Pitch Product of Inertia in Stability Axes 
Roll-Yaw Product of Inertia in Body Axes 
Roll-Yaw Product of Inertia in Stability Axes 


vii i 




M 


M.A.C. 




m. 


N 


q 

r 

[s x ] 
[S 2 3 
T 



Inertia Matrix in Body Axes 
Inertia Matrix in Stability Axes 
Component of [S^] Matrix 
Wing Lift (Force) 

Rolling Moment (Stability Axes) 

Section Rolling Moment (Stability Axes) 

Section Lift Force 

Pitching Moment (Stability Axes) 

Mean Aerodynamic Chord 

Section Pitching Moment (Stability Axes) 

Section Lift Slope 

Yawing Moment (Stability Axes) 

Section Yawing Moment (Stability Axes) 

Roll Rate 

Pitch Rate, Dynamic Pressure 
Rate of Yaw 

Aerodynamic Influence Matrix 
Structural Deflection Matrix 
External Moments 

x-Axis Velocity Component (Stability Axes) 
Flight Air Speed 
Cruise Airspeed 
Downwash Velocity 

Downwash Angle at Three-Quarter Chord Point 
z-Axis Velocity Component (Stability Axes) 
Section Horizontal Force (Stability Axes) 


ix 



X 

Y 

y 


y MAC 

z 


Moment Arm about Longitudinal Stability Axis 
Section Side Force (Stability Axes) 

Moment Arm about Lateral Stability Axis 
Moment Arm about Quarter-Chord Line at A = 0 
Section Vertical Force (Stability Axes) 


C£ Perturbation Angle of Attack 

QIq Cruise Angle of Attack 

a Section Total Angle of Attack 

p Sideslip Angle 

y 0 (Figure 3.4) 

F Strength of Line Vortex 

A Increment (except whenAe(y)) 

Ap Cruise Yaw Angle of Wing 

< 

Sg(y) Flat Wing Section Downwash Angle [rad] 

Ae(y) Section Downwash Angle due to Twist [rad] 

0 Pitch Attitude Angle 

9(y) Section Twist 

cp(y) Wing Built-in Dihedral (unyawed wing) 

P Air Density 

Natural Frequency of Airplane Response Mode 


x 



LIST OF FIGURES 


Figure Page 

2.1 Total Lift and its Component 7 

2.2 Body Axes 10 

2.3 Stability Axes 10 

2.4 Vortex Pattern Representing a Lifting Wing .... 20 

2.5 Horseshoe Vortex Pattern (unya.wed wing) 22 

2.6 Bound Vortex normal to Flight Stream ....... 22 

2.7 Bound Vortex parallel to Wing Span 22 

2.8 Final Angle of Attack 27 

2,9a. Spanwise Lift Distribution 31 

2.9b " " " 32 

2.9c " " " . 33 

3.1 Increment in Load Distribution due to a Change in CL 38 

3.2 Side Force due to Induced Drag 40 

3.3 Side Force due to Wing Rotation 45 

3.4 Effects of Flexibility and Dihedral on Side Force . 45 

3.5 Geometry of Wing Rotation and Moment Arms ..... 51 

3.6 Influence of Wing Rotation on Ailerons Geometry » . 52 

3.7 Perturbation Angles 55 

3.8 Section Lift and Drag Force for an 0! Perturbation 56 

3.9 Sideslip 60 

5.1 NASA-BOEING Design for Oblique-Winged Transport . . 77 

5.2 Total Lift 79 

5.3 Downwash Angles . 81 


xi 



Page 

5.4 Side Force 82 

5.5 Flat Wing Induced Drag 82 

5.6 Ax due to a. Perturbation 84 

5 . 7 Ay due to Q! Perturbation 84 

5.8 Az due to Rate of Roll 85 

5.9 Az due to Rate of Pitch 86 

5.10 Total Lift for Sideslip 87 

5.11 Total Lift 90 

5.12 Flat Wing Lift 91 

5.13a Downwash Angle 93 

5 . 13b Flat Wing Induced Drag 93 

5.14 Total Lift for CL- Perturbation 95 

5.15a Ax due to a-Perturbation 97 

5.15b Ay due to OPerturbation 97 

5.16 Twist 99 

5.17 Total Lift for Roll Perturbation 1°° 

5.18 Az for Roll Perturbation 1°1 

5.19a Ax for Roll Perturbation 1° 8 

5.19b Ay for Roll Perturbation 1°3 

5.20 Total Lift for Pitch Perturbation . 104 

5.21 Az for Pitch Perturbation ....... 105 

5.22a Ax for Pitch Perturbation 107 

5.22b Ay for Pitch Perturbation 107 

5.23 Total Lift Distribution for Negative Sideslip Simulation 1°8 

5.24 AZ for Negative Sideslip 109 

xii 



Page 


5.25a, AX for Negative Sideslip 

5.25b Ay for Negative Sideslip m 

5.26 Total Lift for only Left Aileron Deflected ..... 115 

5.27 Total Lift for only Right Aileron Deflected .... 116 

5.28 Total Lift for Antisymmetric Aileron Deflection . . 117 

5.29a Az for Left Aileron Deflection only . . 118 

5.29b AZ for Right Aileron Deflection only 118 

5.29c AZ for Antisymmetric Deflection 118 

5.30a AX for Left Aileron Deflection only ........ 119 

5.30b AX for Right Aileron Deflection only 119 

5.30c AX for Antisymmetric Deflection 119 

5.31 Total Lift for Left Aileron Deflection only 

(E.A. not Coinciding with Quarter-Chord) 121 

5.32 Total Lift for Right Aileron Deflection only 

(E.A. not Coinciding with Quarter- Chord) 122 

5.33 Total Lift for Antisymmetric Aileron Deflection 

(E.A. not Coinciding with Quarter- Chord) 123 

5.34 AZ for Aileron Deflection 124 

5.35 Lift coefficient, lift curve slope, and 

aerodynamic center versus skew angle 143 

5.36 OL stability derivatives versus skew angle 144 

5.37 /3 stability derivatives versus skew angle 145 

5.38 p stability derivatives versus skew angle ...... 146 

5.39 q stability derivatives versus skew angle 147 

5.40 r stability derivatives versus skew angle 148 

5.41 u stability derivatives versus skew angle 149 


xiii 



5.42 Root locus versus skew angle for 

spiral mode 

5.43 Root locus versus skew angle for 

rolling convergence 

5.44 Root locus versus skew angle for 

phugoid mode 

5.45 Root locus versus skew angle for 

short period mode 

5.46 Root locus versus skew angle for lateral 

oscillation or Dutch roll . . . 

5.47 Effect of skew on rolling convergence 

and spiral mode shapes 

5.48 Effect of skew on the phygoid mode shape 

5.49 Effect of skew on short period mode shape 

5.50 Effect of skew on lateral oscillation on 

Dutch roll 

5.51 Rigid Wing 40$ SAS - Response to Aileron 

5.52 Elastic Wing 40$ SAS - Response to Elevator . . . . 

B.l Sideslip of an Oblique Wing Aircraft 

B.2 Moment about the C.G. in the Plane of Symmetry. . . 

D. 1 Horseshoe Vortex System 

D. 2 Downwash Control Points 

D.3 Finite Segment of a Straight Vortex Filament. . . . 

D. 4a Control Point to the Right of Horseshoe Vortex . . 

D.4b Control Point Within the Horseshoe Vortex . . . . 

D. 4c Control Point to the Left of Horseshoe Vortex . . . 

D. 5 Structural Skeleton of the Wing . . . 

D.6 Plan View of Left Wing Tip Section 

D.7 Aerodynamic Moments Balance at the Root 


Page 

153 


153 

154 


155 


156 


157 

158 

159 


160 

162 

163 

171 

176 

191 

193 

196 

198 

200 

202 

206 

209 

211 


xiv 



D.8 Plan View of Aircraft Centerline . 214 

D. 9 Structural Moments due to Control 226 

D. 10 Plan View of the Geometric Parameters used in 

Computing the Apparent Twist due to Angular 

Velocities . 230 

D. 11 Wing Section Center of Mass 234 


xv 



LIST OF TABLES 

Page 

Laplace Transform of the Equation of Motion .... 72 

Ailerons Dimensions 113 

Influence of the Stability Derivatives on the 

Natural Modes 128 

Comparison Between Stability Derivatives , . . . . 130 

Computer Printout for Case 1 of Table 5.3 ........ 132 

Computer Printout for Case 2 of Table 5.3 133 

Computer Printout for Case 3 of Table 5.3 134 

Computer Printout for Case 4 of Table 5.3 ..... 135 

Computer Printout for Case 5 of Table 5.3 ..... J36 



I - INTRODUCTION 


1. 1 Background , 

The oblique wing solution for a transonic transport has been studied 
increasingly in the past few years. This is not a new idea; it was 
first proposed more than 25 years ago and, in addition to work in this 
country, it appears that German aerodynamicists made classified studies 
of this configuration during World War II, although the project never 
reached- the stage of flight tests [Ref. 1]. 

In 1947, a first attempt to investigate the stability of an oblique 
wing was made in the NACA Free Flight Tunnel. The results of the in- 
vestigation indicated that is was possible to skew the wing as a unit 
to angles as great as 40° without encountering serious stability and 
control difficulties. At an angle of skew of 60°, however, the aile- 
ron control became unsatisfactorily weak [Ref. 2]. 

For supersonic speeds, in addition to the friction and vortex 
drag, the wing experiences a wave drag associated with the thickness 
or volume of the wing, as well as with the lift distribution. The 
vortex drag is independent of the distribution of lift in the flight 
direction, whereas the wave drag does • depend on the distribution in 
the direction of flight. However, it diminishes as the length of the 
wing in the flight direction increases. In fact, linear theory shows 
that the wave drag due to lift diminishes approximately as the inverse 
square of the length, while wave drag due to volume diminishes as the 
inverse of the fourth power; furthermore, spreading the lift over a 
greater length diminishes the sonic boom intensity [Ref. 3]. 


1 



Structurally, an oblique wing presents the advantage of having a 
structure continuous across the pivot and makes tension the primary 
load on the pivot. It is well known that swept- forward wings show a 
tendency for aeroelastic divergence and this remains an area of concern 
for the oblique wing also. In any case, recent Boeing studies show that 
for the unrestrained airplane it is possible to reach a static stabil- 
ity without a severe penalty in weight. Oscillatory aeroelastic in- 
stability occurs at speeds higher than those at which the clamped 
fuselage static instability occurs (the clamped fuselage being a more 
conservative approach) [Ref. 4]. 

The object of this study is to analyze the aerodynamics and 
dynamics of oblique wing aircraft. 

So far, with the exception of limited work based on extensions of 
symmetric aircraft theory to the yawed wing aircraft concept, only 
Boeing Commercial Airplane Company and NASA-Ames Research Center have 
done systematic studies on the dynamics of oblique wing aircraft. The 
study by Boeing [Ref. 25], a NASA sponsored one, was aimed toward a 
design investigation for flight in high transonic speed regime. The 
yawed wing aircraft was chosen as the most promising configuration for 
that regime. During this work, we shall use the design solution for 
a yawed wing aircraft chosen by Boeing and compare portions of our 
results with theirs. 

The work done at NASA-Ames Research Center in the last few years 
has been characterized by several wind tunnel tests to investigate wing 
aerodynamic coefficients at different skew angles. Future programs are 
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very interesting. A remote-pilot vehicle (RPV) with an oblique wing, 
having 26' (unyawed) span, has already been built. The flight tests, 
investigating low subsonic speed, are expected to start before the end 
of the year. A pressure model, based on Boeing design 5-3 [Ref. 25] 
should be completed within two years. 

A supersonic drone, the Firebee II, a Flight Research Center (FRC) 
project capable of Mach 1.4, should be completed in three years. The 
data gathered from this project will be used for the one which in- 
volves the modification of a NASA fighter, the F-8, with an oblique 
wing. This last project, another FRC program, should be completed 
within three years. 

Very recently, Boeing has considered an oblique wing design for 
a commercial transport cruising at Mach .95. 

1.2 Synopsis and Contributions . 

As an mid to the reader, the remaining chapters of this report 
are briefly summarized below. 

Chapter II deals with the wing aerodynamics. Spanwise distribution 
of lift and induced drag are computed according to the strip theory and 
lifting line theory methods. An empirical correction to Schrenk's 
method for evaluating spanwise lift distribution [Ref. 6] is then 
proposed for the oblique wing. 

In Chapter III, the basic features of an oblique wing aircraft are 
first analj^zed. A lift contribution to the side force due to wing 
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rotation has been found and its derivation is described; the induced 
drag also introduces a side component and its magnitude is computed. 

The new derivatives due to the oblique wing and the general methodol- 
ogy for evaluating stability derivatives are also described in this 
Chapter. 

Chapter IV deals with the derivation of the equation of motion. 

The numerical results are given in Chapter V. First, the span- 
wise distribution of lift, downwash, and induced drag, as computed by 
strip and lifting line theory, are reported and compared. Then, the 
results for rigid and flexible wing, as computed by lifting line theory, 

are described. Special emphasis is given to the spanwise distribution 
of the increments of the aerodynamic forces during simulated perturbed 
conditions, and the resulting stability derivatives. Root loci vs. 
skew angle along with the mode shapes are shown for all the natural 
modes. 

Chapter VI contains the final conclusions and remarks. 
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II AERODYNAMICS 


2. 1 Generalities 

The lift generated by a wing can be divided into two contributions: 
lift due to twist and camber, and lift due to the flat wing at angle of attack. 

According to linear theory these two contributions are independent 
of each other* and the lift is simply equal to their superposition. In 
addition, the twist contribution does not depend on the angle of attack 
and therefore on its variations, but it does depend on sweep and speed 
changes (or, more exactly, on dynamic pressure variations). In a 
stability analysis we must know the spanwise distribution of the aero- 
dynamic force and how such distribution varies when the wing undergoes 
perturbations. Let us now limit ourselves to the lift component of the 
aerodynamic force since drag and side forces can be derived from the 
knowledge of the lift. 

For a symmetric aircraft it is possible to obtain the lift distri- 
bution for the perturbed condition by applying strip theory to the cruise 
spanwise lift distribution [Ref. 7 ]. Though an approximate one, this 
method produces satisfactory results in the linear range. 

The perturbations considered in a stability analysis can be divided 
into two groups: 

1) perturbations which affect the total lift distribution (flat wing 
plus twist contribution) 

2) perturbations which affect only the flat wing. 

To the first group belong the u , r , and |3 perturbations; to 
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the second group the perturbations involving change in the local angle 
of attack (a , p , and q) . 

For a symmetric aircraft the flat wing spanwise lift distribution 
for the cruise condition is symmetric. 

Let us now consider the flat wing lift distribution (FWLD) of a 
skewed wing. For such a wing at no sweep, the lift distribution is 
symmetric and elliptic (because of the chord distribution). When the 
wing is swept, the new spanwise lift distribution will now show an 
increase in lift on the aft part of the wing and a decrease on the 
forward one, because of an upwash generation in the first case and a 
downwash increase in the second one. The need for twist in this wing 
will be based now not on a ’tip stall last’ requirement, but primarily 
on the need to obtain a symmetric lift distribution from the asymmetric 
flat wing lift. Fig. 2.1 shows the asymmetric lift, lift due to twist, 
and total lift versus span. The accurate evaluation of the change in 
aerodynamic loading will lead to good precision in computing the wing 
contribution to stability derivatives. We shall now follow two approaches 
in determining the spanwise lift distribution during the perturbed motion. 
The first, based on strip theory, can be used for all perturbations, but 
in the case of sideslip it fails to give any acceptable result. In 
order to improve this situation an empirical method, based on a modifica- 
tion of the Pope-Schrenk method [Ref. 6] which can compute the spanwise 
lift distribution for a flat oblique wing was obtained and will be 
briefly discussed in section 2.5 , 

The second approach is based on the evaluation of the spanwise lift 
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distribution of an oblique wing by means of linear aerodynamic theory. 

In this case, the spanwise lift distribution is computed for both cruise 
and perturbed conditions. This second method also allows the static 
effects of a flexible wing to be included. The load distributions ob- 
tained by means of these two methods will then be used in determining 
the stability derivatives and the results compared. Once we have com- 
puted the stability derivatives we shall use these values in the equations 
of motion and determine the natural behavior and the time response to 
control disturbances. 

2.2 Reference Axes 

To discuss the problem of stability, it is necessary to set up a 
system of reference axes which form the basis of a system of notation 
used to describe the motions of an airplane. Two basic body fixed axes 
systems, each consisting of three mutually perpendicular axes passing 
through the center of gravity of the airplane, adequately cover most of 
the aerodynamic problems in stability considerations. These are the 
"body axes" and the "wind axes", which will be referred to as "stability 
axes". 

2.2.1 Body Axes (x^ , y^ , z^). 

The body axes system is rigidly fixed in the airplane and is the 
system of mutually perpendicular axes passing through the airplane's 
center of gravity and whose X-axis is parallel to the thrust axis, the 
wing mean aerodynamic chord, or some other longitudinal reference and is 
positive in the direction of the nose of the airplane. Herein it is 
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taken along the body centerline reference. Figure 2.2 shows this 
system together with all forces, moments, displacements, and velocities. 
The XZ plane is the plane of symmetry for the airplane. 

2.2,2 Stability Axes (x^ , y^ , z^). 

The stability axes system differs from the body axes system in that 
the X-axis is parallel to the relative wind, positive forward. The 
Z-axis is positive down (Fig. 2.3). The Y-axis is positive to the right. 
The moment, angle, and angular velocity conventions are given by the 
right hand rule. 

The stability axes system is the one used for basic aerodynamic 
performance work, and it will be used in our analysis. In our study we 
shall also assume that for zero angle of attack stability and body axes 
will coincide. 

2 . 3 Strip Theory 

The strip hypothesis asserts that we may calculate the aerodynamic 
force on each strip as if it were an isolated airfoil moving with the 
resultant velocity which it has because of its local position on the 
aircraft [Ref. 7]. 

The strip hypothesis is a first approximation to the actual case 
where the trailing vortices from each strip possibly interfere with the 
others, but this first approximation is confirmed by experiments to give 
results which are in excellent agreement with the facts as observed, 
even at incidences above the stall, for a. symmetric aircraft. Unfortu- 
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Notation for body axos. 


L — rolling moment p — rate of roll 

M = pitching moment q = rate of pitch 

N == yawing moment r - = rate of yaw 

[X, Y, Z] — components of resultant aerodynamic forco 
[u, v, w] = components of volocity of C rolativo to atmosphere 

Figure 2.2 - Body Axes [Ref. 21] 
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nately, when dealing with motions involving rate of pitch or rate of 

roll, the strip theory fails to give a good approximation for the case 

of an oblique wing, since it predicts no lift change at the wing root 

chord for the case of a roll, or at the wing station, whose quarter 

chord is crossing the y axis, for the case of pitch, 

s 

The advantage of this approach is that it allows us to evaluate 
the stability derivatives using different sources of data, such as 
wind tunnel or flight tests, and computer results, 

2.3.1 Variations in Aerodynamic Forces due to a Perturbation a . 

Since the twist (or dihedral) contributions are independent of 
changes in angle of attack, only the flat wing lift distribution (FWLD) 
needs to be considered. (Subscript F will indicate flat wing quantity). 

We can now introduce the following approximations. 

2.3. 1.1 Section Lift Slope. 

Assuming the FWLD is given, the local lift slope can be computed as 


where 


c T (y) 


[C L (y)] F 


( 2 . 1 ) 


C (y) ~ local lift slope at station y (1/ra.d) 

L a 

[C T (y) ] = local lift coefficient of the flat wing at 

L r 

station y 


Qq = angle of attack for cruising condition (rad). 
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2.3. 1.2 Section Downwash. 

The downwash is a consequence of the wing not having an infinite 
span. For a 2-D wing, the section lift coefficient is given by 

t C L (y) 2-D> Tot = ra 0 [a 0 + 6(y)] = m 0 a 0 + m 0 6(y) (2 ' 2) 

where 

mg = 2-D section lift slope 
0(y) = section twist. 

The presence of trailing vortices in a finite wing introduces local 
downwash velocities whose effect is to reduce the local angle of attack 
and, therefore, the lift produced by the wing. The two terms of (Eq . 
2.2) will become, for a 3-D wing 


a o“o = »o[“o'V y)] 

a Q 0(y) = m 0 [9(y) - A«Ky)] 

where 

Cg(y) = flat wing downwash angle 

Af(y) = downwash angle due to twist. 


(2.3) 


Therefore, the section total lift coefficient for a 3-D wing is 


CC L (y)] Tot = m oj [ V V y)] + te(y) - A€(y)]J (2.4) 

The previous expression applies to a straight as well as to a skewed 
or swept back wing; the only change, assuming all quantities are measured 
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in the flight direction, would occur in the value of a^ , the 2-D 
section lift slope. We refer to R.T. Jones [Ref. 7] for a detailed 
discussion of the derivation of a.^ . 

We shall now evaluate the downwash angles by computing the difference 
in the local lift distribution from the 2-D case. This is a "crude" 
approximation, but in a strip theory analysis it is the only way to 
evaluate the downwash velocity and, therefore, the spanwise induced 
drag distribution. 

In section 5 . 2. 1 we sha 11 compare the downwash results obtained with 
this method against the corresponding results obtained applying the 
method based on linear theory described in the next section. 

Let us define 

t C L (y)] Tot= t C L ( ^F + (2 ' 5) 

where 

[C T (y)]_ = flat wing section lift coefficient 
L F 

[C_ (y) ]_ = section basic-lift coefficient. It represents the 
L B 

a- independent contribution to lift due to twist, 

[c L (y)] F = a 0 [« 0 " V y >] ( 2 * 6 ) 

[C L ( y) ]B = V e < y) - Ae ^ (2.7) 

We can now assume, according to linear theory, that the two down- 
washes are independent from each other, and furthermore, that the down- 
wash due to twist does not depend on the angle of attack. With these 
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assumptions we have 


[, C t (y) ]f = a n[ a n - = a n ) 1 “ 


9e Q (y) 


'L w/j F j 0 l ~0 ' 0 w/j “0 )~ a . a o 


( 2 . 8 ) 


d€ 0 (y) _ ^ _ [ c L <y)] F _ i _ l \ (y)1 


da 


a o a o 


0 


(2.9) 


de Q (y) 

*0 (y) = ~da~ a Q 

The total downwash is therefore given by 


( 2 . 10 ) 


[e( y ) ] Tot = € Q (y) + Ae(y) 


( 2 . 11 ) 


9Le( y ) ] Tot ^ 9^ 0 (y> A^( y ) 1 L a 

3 a 3a 3a a Q 


( 2 . 12 ) 


2.3. 1.3 Section Induced Drag. 

The aerodynamic drag in a finite wing has two components: the first 
is due to skin friction and pressure distributions on the boundary, the 
second is the one induced by the lift because of the presence of trailing 
vortices. 

Both components are normally of the same order, and dependent on 
the aircraft speed [Ref. 15]. The analytic spanwise evaluation of these 
two components is a difficult task. 

According to Multhopp [Ref. 9] we may write the induced drag coef- 
ficient for the wing as 
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b/2 

Id . 2s / 

1 -b/2 


c C i/ y ^ a i dy 


(2.13) 


where the so-called induced incidence is 


b/2 

It 

-b/2 


_d_ 
?7 dTj 


CC L (y) 


dt] 


(2.14) 


and the section induced drag coefficient is given by 


C D (y) = C L (y)a i (2.15) 

1 

Garner, [Ref. 10] who has discussed induced drag and its spanwise 
distribution in incompressible flow, has shown that, for swept back 
wings, the quantity C C_(y) (X, does not give an acceptable spanwise 
distribution of induced drag as suggested by Robinson and Laurmann 
[Ref. 11], This conclusion can be expected since the induced incidence 
Ct , as computed in Eqn. (2.14) implies that all the bound vortices of 
the horseshoe vortex system, used as a model for the wing, lie on a 
straight line perpendicular to the velocity. Therefore, the downwash 
angles correspond to the ones of a straight wing and the induced drag 
obtained is the product of the lift distribution of a swept wing times 
the downwash angle of a straight wing having the same wing span and the 
same spanwise lift distribution. 


We shall return to this subject when we evaluate the induced drag 
distributions by means of linear theory. For the sake of clarity we 


w 


shall recall that the downwash angle is the ratio — where 

V 
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w = downwash velocity measured a.t the lifting line 
V = free stream velocity. 

In our strip analysis the downwash angles are obtained directly 
from the actual spanwise lift distribution; we assume that the downwash 
angle is the cause of the difference in lift coefficient from the 
corresponding 2-D one, and no assumptions are made on the wing geometry. 
We may therefore expect a. better accuracy in the estimate of the spanwise 
drag distribution. 

Let us now see how to relate the previous discussion to the evalua- 
tion of stability derivatives by means of strip theory. 

For the lift case we have one term, the basic lift, which is 
"a independent" and the other, the flat wing, which does depend on cc ; 
the same considerations can be applied to the drag. In the "Q! indepen- 
dent" drag contribution we can group skin friction, pressure distorsions, 
and, if we extend to the drag the same assumptions made for the lift, 
also the induced drag due to basic lift. This approximation is quite 
accurate for skin friction and pressure distortion, but it becomes less 
accurate when considering the induced drag of the basic lift, as is 
shown next. 

The total section induced drag is obtained by substituting into 
Eqn.(2.15) the total values derived in the previous sections 

[ c D _(y)] To t “ ^ C L < - y ' ) ^Tot^ C 0^ \ot 
1 


16 



(2. 16) 


- |[C L (y)] F + [c L (y)] B j t€ 0 (y) + Ac(y) ] 

= [C L (y)l F [c Q (y) + &(y)] + [C L (y)] B [€ Q (y) + Ac (y) ] 

We can see that the induced drag produced by the basic lift 
distribution 

[C D (y)] = [C L (y)] B [€ 0 (y) + Ac (y) ] (2.17) 

i 

has the term [C (y)] € n ( y) that is a dependent since C n (y) , as 

Li B U u 

discussed in the previous section, is (X dependent. 

Therefore, only the quantity [C (y)] AC (y) can be assumed, 

L B 

within the range of linear theory, as a independent, whereas 

[C L (y)] B V y) + tC L (y)] P [€ 0 (y) + AC(y)] (2 ' 18) 


is the a dependent component of the drag. The change in the section 
drag for perturbations introducing a local change in the angle of attack 
can now be computed as the rate of change of Eqn. (2.18) with CX 


3C d (y) 

-ir- ■ i e o« + t C L<^F + 


[C T (y)]. 


3e 0 ( y ) 


(2.19) 


+ C T (y) [€ n (y) + Ac (y) ] + C T (y) ], 
‘a 


and by substituting Eqns. (2.5), (2.11) and (2.12) into (2.19) we obtain 


SC D _(y) 

to ~ ^L^^Tot 


1 - 


c T (y) 

L a 


a 


0 1 


+ C L (y) [e 0 (y)] Tot (2.20) 
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2.3.2 Evaluation of Stability Derivatives by Means of Strip Theory. 

Assuming the spanwise lift distribution is given for the flat wing 
as well as for the wing with the nominal twist and dihedral, it is pos- 
sible to evaluate the stability derivatives by using the approximations 
used in section 2.3.1 and the expressions given in sections 3.3.2 and 
3.3.4 . 

For the stability derivatives due to side-slip it is necessary to 
know the spanwise lift distribution corresponding to the new sweep angle 
differing from the nominal one by A = - p . Because of the peculiarity 
of the shape of the FWLD and its dependence on the sweep angle, an ap- 
proximation assuming 

( [CC LW\o t )A- (2 - 21) 

13 cos Aj, 

where 

A = A - p 

would fail to give an acceptable result. The knowledge of the spanwise 
lift distribution for the new sweep angle would therefore be required. 

For the case when such lift would not be available, an empirical 
correction to the Pope-Schrenk' s method [Ref. 6,12] was derived and it 
is described in section 2.5 . 
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2.4 Lifting Line Theory 


2.4.1 Introduction 

A more accurate way of evaluating the stability derivatives, spe- 
cially roll and pitch derivatives, is to use linear aerodynamic theory. 

The literature offers a wide variety of methods which can be used; many 
of these are very complex and allow the user to evaluate both the span- 
wise and chord-wise load distribution. The simplest three-dimensional 
wing theory is that based on the concept of the lifting line. In this 
theory the wing is replaced by a straight line [Ref. 13]. The circula- 
tion about the wing associated with the lift is replaced by a vortex 
filament. This vortex filament lies along the straight line; and at 
each spanwise station, the strength of the vortex is proportional to the 
local intensity of the lift. According to Helmholtz's theorem, a vortex 
filament cannot terminate in the fluid. The variation of vortex strength 
along the straight line is therefore assumed to result from superposition 
of a number of horseshoe-shaped vortices, as shown in Figure 2.4 , The 
portions of the vortices lying along the span are called the "bound 
vortices". The portions of the vortices extending downstream indefinitely 
are called the "trailing vortices". 

The effect of trailing vortices corresponding to a positive lift is 
to induce a downward component of velocity at and behind the wing. This 
downward component is called the "downwash". The magnitude of the down- 
wash at any section along the span is equal to the sum of the effects 
of all the trailing vortices along the entire span. The effect of the 
downwash is to change the relative direction of the air stream over the 
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section. 



Figure 2.4 - Vortex Pattern Representing a Lifting Wing 

[Ref. 13] 

The section is assumed to have the same aerodynamic characteristics 
with respect to the rotated air stream as it has in normal two-dimen- 
sional flow. The rotation of the flow effectively reduces the angle 
of attack. Inasmuch as the downwash is proportional to the lift coef- 
ficient, the effect of the trailing vortices is to reduce the slope of 
the lift curve. The rotation of the flow also causes a corresponding 
rotation of the lift vector to produce a drag component in the direction 
of motion. 

The methods using discrete vortices to represent the continuous 
distribution of circulation of the vortex sheet are attempts to simplify 
the computations. In the methods employing discrete vortices, two-dimen- 
sional theory is used to determine the most representative locations of 
the vortices as well as of the control points. If only one vortex line 
is used, it is placed along the center-of-pressure line in two-dimensional 
flows, which for a flat plate at an angle of attack and at subsonic speeds 
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is the quarter chord line. 


For subsonic speeds, the downwash varies inversely with the 
distance behind the quarter chord line; at the position of the three- 
quarter chord line it just equals in magnitude the vertical component 
of the flow tangential to the flat plate having the same circulation. 
Conversely, if the condition of tangential flow is satisfied at the 
three-quarter chord line, the strength of the concentrated vortex will 
indicate the lift on one wing due to angle of attack. Of course, these 
methods of obtaining correspondence between the lifting lines and the 
vortex sheets loose their validity near the corners of the wing, where 
the flow differs sharply from the two-dimensional. 

In a stability analysis it is not required to have great accuracy 
in the chordwise load distribution. Therefore, it is sufficient to use 
only one vortex line placed, for subsonic cruise condition, at the 
quarter-chord point. In doing so, little is lost in accuracy and a lot 
is gained in simplicity. 

The supersonic case has to be approached in a different way accord- 
ing to the properties of the supersonic flow. The method outlined in 
this chapter applies to subsonic speeds only, but it can be extended 
to the supersonic case. It can be shown that the horseshoe-vortex 
system of Figure 2.4 is equivalent to the one where each horseshoe vortex 
has a constant strength equal to the sum of the strength of the bound 
vortices contained at the corresponding section of Figure 2.5 . 

When the wing is skewed, two models can be used. Figure 2,6 shows 
the first one where the bound vortices are aligned with the wing span , 
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and are at a skew angle with, respect to the free stream velocity. 

The second model, shown in Figure 2.7, assumes that the bound 
vortices are perpendicular to the flight direction. This model corre- 
sponds to the case where the wing at a skew angle is actually replaced 
by a finite number of straight wings, and will be used in this study. 

The first model is expected to give a better accuracy when the 
number of horseshoe vortices is small, but the two models coincide when 
the number of horseshoe vortices goes to infinity. 

The method used in this study is based on a modification of the 
Weissinger-L-Method and applies at subcritical Mach Number. This 
method is derived from excellent work done by Gray and Schenk in 1953 
[Ref. 14] and has been modified for the oblique wing case. Among the 
advantages of such an approach to the evaluation of the spanwise dis- 
tribution of the loading, is the possibility of evaluating the loading 
increments due to aileron or flap deflection, effects of wing flexibility, 
and accelerations on the wing. These features are fully described in 
Appendix D. 

2.4.2 Steady State Loading on an Airplane with an Oblique Wing. 

The fundamental problem involved is the development of a series of 
equations which relate the spanwise lift distribution for an arbitrary 
wing plan form in a given flight condition to the properties and atti- 
tudes of the individual sections that form the wing. 

For a 2-D wing, the following relationships can be found in any 
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standard textbook on aerodynamics [Ref . 7 , 15 » 17 ] : 


where 


i = pvr = q tn 0 C a 



m^ = local lift slope 
C = local chord length 
i = section lift 

C£ f - total angle of attack (see Fig. 2.8) 
1 „ 2 , 

q = — p V = dynamic pressure 
r = circulation 


( 2 . 22 ) 

(2.23) 


At a specific distance r behind the lifting line, the resultant 
of the downwash velocity a.nd the flight velocity V is parallel to 

the section zero lift line. Then, 


W = CL. V 
r f 


(2.24) 


from Eqn. (2.22) we obtain 


m 0 q C a 

r - -V 4 - “0 f C “f 


(2.25) 


Substituting (2.25) into (2.23) results in 


_ m 0 C/2 

W r - ~ “f V 


(2.26) 


Equating (2.26) and (2.24) 


V - 3} SO. a. V 
f 2TT r f 


or 


m 


_0 C/2 
277 r 


= 1 


(2.27) 
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Since the theoretical section 2-D lift curve slope is 2ff , r must 


equal C/2 , which is the distance between the lifting line and the 
three-quarter-chord point. 

Therefore, for the 2-D (unswept wing) ca.se 

(1)30/4 * “ f < 2 - 28) 

corresponds to the control point where no flow exists normal to the zero 
lift line. Whenever the local lift slope differs from 2 tt , expression 
(2.28) becomes 

(l) 3 c/4 * “f (2 - 29 > 


The essential difference between a 2-D wing and a wing of finite 
aspect ratio arises from non uniform spanwise loading which produces 
the trailing vortices. The equations presented so far are considered to 
apply to a finite wing when the effects of all the vortices, both bound 
and trailing, have been taken into account. 


Equation (2.29) can be written in matrix form 


|W( 

(v(3C/4 




(2.30) 


This matrix relation represents a series of equations, each appli- 
cable to a. particular station on the span of the wing. 

The elements of | v ! 3 /A C ’ ever ^ one which is affected by the 
bound and trailing vortices of the wing stations can be evaluated from 


(W 

I V )3’/4 c 


1 

4 7TV 


[SI] r 


(2.31) 
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and by expressing P in terms of \i 


The [Si] matrix is the downwash matrix and is derived in Appen- 
dix D. 


Combining equations (2.30) we obtain 









or 


1 

4qm 0 


[SI] 




(2.32) 


(2.33) 


2.4.3 Section Final Angle of Attack { | . 

The final angle of attack across the span jcikj can be considered 
to be composed of three essential parts (see Fig. 2.8). 






where 


(2.34) 





= angle of attack caused by structural deflection of a 
flexible wing 

= angle of attack caused by built-in twist, apparent or 

aerodynamic twists, control deflection, angular velocities. 
= flat (and rigid) wing angle of attack (measured w.r.t. 
root-section zero-lift line). 


The angle of attack Q! g j-, caused by structural deflection of a 
flexible wing due to the section lift at the section aerodynamic centers 
is linearly related to the matrix jij as 
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Figure 2.8 - Final Angle of Attack 


op's'll 
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(2.35) 


! a ,! - ts 2 ] Ml 

where [S^] = structural deflection matrix (described in Appendix D). 

The contributions to the angle of attack Ct are described in 

g 

Appendix D. 


2,4.4 Section Induced Drag. 

As discussed in section 2.3,1, only the induced drag varies with a. 
The induced drag arises from the rotation of the aerodynamic forces due 
to the downwash velocities induced by the trailing vortices. 


In our model, where the wing has been replaced by a lifting line 
placed, for the subsonic analysis, at the quarter-chord line, the 
section induced drag can be evaluated by computing the downwash angle 
(OjJi at the station bound vortex. This is similar to what is done 


for * 


J/3/4 C 
is zero. 


, except that downwash induced by a bound vortex on itself 


Therefore 



TO 


E i ( Vc/4 i) - 


PV, 


TO 


^j (K ij } C/4 


i. 

J 


(2.36) 


where (Kij) c ^ is the same as computed in Appendix D with the following 
two exceptions: 

C. 

1 

1) the terra — , distance of the control point from the lifting 
line, must be dropped since the downwash is now evaluated at 
the lifting line; 
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2) the contribution to the downwash sould be disregarded for the 
case when the control point is within the horseshoe vortex, as 
discussed in Appendix D, 

A more accurate result would be obtained using the vortex lattice 
method. The results obtained in this analysis are discussed in Refer- 
ence 16 . 

2.5 An Empirical Correction to Schrenk's Method 

Another approach to the problem of spanwise load distribution -t 

having much less theoretical foundation is presented by Flatt [Ref, 18]. 

It follows a method first presented by Schrenk [Ref. 19]. 

Schrenk's method makes allowance for the effect of the varying down- 
wash along the span of a nonelliptic wing by assuming that the final 
span load distribution for an untwisted wing is halfway between the 
actual planform shape and a semi-ellipse of the same area. However, 

Schrenk and Flatt did not consider a swept back wing; Alan Pope and 
William R. Haney, assuming that the effect of sweepback on the non- 
dimensional span loading is linear, proposed an empirical correction 
[Ref. 6, 12] to take care of the effect of sweep back. It has been 
successfully employed in preliminary design of subsonic aircraft. 

The following empirical formula which was obtained during our 
study of the oblique wing, can be applied to Schrenk's method. 



29 



where 


A = skew angle 

c =■ c /cosA : root chord measured in the flight direction 
R R q 

b = cosA : actual wing span 

b Q = wing span for the unskewed wing. 

Since the goal of this method is to provide, in a fast way, the 
shape of the spanwise lift distribution, the wing lift coefficient is 
assumed to be known. In addition, the new spanwise lift distribution 
obtained with this correction must be normalized by multiplying it times 
the ratio of the wing lift coefficient to the lift coefficient obtained 
by integration of the new lift distribution. 

The results of this correction, for the case of a flat wing, were 
checked against the lift distribution as computed by the numerical 
program of Reference 20. 

Figures 2.9a, b, and c show the comparison between the two methods 
for three different cases. 

To date, no check has been done for the case of a. wing having 
twist. The major difficulty for this case is to estimate the 
downwash velocities produced by the lift due to twist and, therefore, 
to determine the effective twist. 

Though for symmetric wings experience suggests that an effec- 
tiveness of 50^ is an acceptable assumption, we do not expect that the 
same can be applied to a shewed wing. 

As pointed out before, the effect of a. positive sideslip corresponds, 
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for an oblique wing, to a negative change in sweep and the perturbation 
will therefore affect both the FWLD and the basic lift distribution. 


To a first order approximation it is possible to assume that the 
change in the shape of the lift distribution for the total lift will 
be proportional to the change with sweep of the FWLD. 

Therefore, the simulation of a sideslip can be done by 

1) determining , with the Pope-Schrenk' s method corrected for the 
skewed case, the FWLD for both the cruise sweep A^ , and the 
sideslip condition A = A^ - [3 ; 

2) multiplying the spanwise distribution of the ratio of the FWLD 
at A . = Aq to the FWLD at A = A^ - p times the total lift 
distribution for the cruise condition in order to obtain the 
corresponding distribution for the sideslip condition. 



2 . 6 Summary . 

The reference axes system used in this work, the stability axes, 
has first been defined. The method for evaluating variations in the 
aerodynamic forces due to perturbations was then outlined. A simplified 
method, based on strip theory, for computing spanwise distribution of 
the induced drag was also proposed. 

A second and more systematic way of computing spanwise lift 
distribution based on lifting line and including effects of wing flexi- 
bility was the described. The spanwise induced drag distribution was 
then evaluated by using lifting line theory. 

In the last section of this Chapter, an empirical correction to 
Schrenk' method for the oblique wing case was proposed. 


35 



Ill STABILITY DERIVATIVES 


3. 1 Basic Features 

3.1.1 Flat Wing Lift Distribution (FWLD) . 

In section 2,1 we have described how the FWLD behaves when the wing 
is skewed. This behavior is undoubtedly one of the most important 
elements of difference from the symmetric case. The upward field gener- 
ated by the forward wing, causes the aft wing to "feel" a higher angle 
of attack. This, as we mentioned before, is the cause of the loss of 
symmetry in FWLD, but, because of the higher angle of attack actually 
experienced by the aft wing, the stall condition will be reached in this 
region first. This is similar to the case of a tip stall for a symmetric 
wing, except that now only one tip would sta.ll and the loss in balance 
would produce not only rolling but also pitching moments. The recovery 
from such a stall is very difficult. 

This can be explained if we analyze in detail the motion of the 
aircraft following the sta.ll of the aft wing. 

As we mentioned before, the loss in the lift symmetry introduces, 
in the case of the left wing forward, positive pitching and rolling 
moments. The nose-up motion deriving from the pitching moment introduces 
an increase in the angle of attack of the whole wing as well as an 
angular velocity; whereas the rolling moment produces only an angular 
velocity. 

The effect of these two angular velocities results in a linear 
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variation of the local angle of attack for the wing, increasing it on 
the aft part and decreasing it on the forward one. 

Therefore, the aft wing will experience a further increase in the 
angle of attack which will worsen the stall condition and extend it 
inboard. 

On the forward wing, instead, the increase in angle of attack of 
the whole wing is counteracted by the decrease deriving from the angular 
velocities and the final trend is consequently toward a delay in the 
forward wing stall, which makes the nose-down attitude, required for a 
recovery, difficult to achieve. 

The center of pressure of the flat wing lift distribution (FWLD) 
lies on the aft- portion of the wing and produces both pitching and 
rolling moments as the angle of attack varies, introducing new important 
derivatives. In fact, any change in the angle of attack affects the 
FWLD only and destroys the symmetry in the total lift distribution 
obtained by twisting the wing. The consequence of this is a loss in the 
moments' balance. Figure 3.1 shows the variation in spanwise lift 
distribution, as given by linear theory, for a case where, for the sake 
of clarity, the "perturbation" a is assumed to be equal to 50$ of . 

It will be seen that a positive change in a produces an increase in 
the lift whose aerodynamic center is displaced toward a point in the 
aft part of the wing, thus producing a rolling moment as well as a 
pitching moment. By considering the corresponding change in induced 
drag, it is possible to evaluate the yawing moment. 
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3.1.2 Side Force 


A wing at a skewed angle experiences a side force deriving from 
the induced drag component in the y direction. Let us now derive 
its magnitude. 

Figure 3.2 shows a section of a wing at a skew angle A . 

The principle of independence [Ref. 8] assumes that, in a friction- 
less flow, all the results of the two-dimensional flow theory can be 
applied immediately to an infinite oblique wing simply by subtracting 
the axial component of velocity. 


No matter whether we consider the free stream direction or the 
direction normal to the leading edge, the lift generated by the wing 
section is obviously the same. Therefore, the lift force is given by 


n 1 „2 „ r ^ 1 „2 „ . C ■ 1 

l ‘i pv o c L 0 pv V^A 


(3.1) 


where the subscript 0 defines the quantity in the flight direction. 


Since 

V = V Q cosA (3.2) 

and 

C = C Q cosA (3.3) 

from 3.1 it is possible to derive the well known relationship 

C = C T cos 2 A (3.4) 

L 0 L 

Figure 3.2a shows the decomposition of the velocity vector into two 
components; one normal and one parallel to the leading edge. The two 
elementary strips have the same area. 
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Figure 3.2 - Side Force due to Induced Drag. 


If we now introduce a plane defined by the velocity vector Vq 
and by the trailing edge segment where the two strips coincide, we can 
derive the relationship between the angles of attack measured in the two 
directions. This relationship is useful when evaluating the angle of 
attack during sideslip. This can be done by simply noticing that the 
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height from the leading edge to the previously defined plane is the 
same at the point as well as at k . This, of course, implies 

that the leading edge is a straight line. The effects of twist or 
dihedral would then be superimposed. 

The distance h being a constant, by pure trigonometric consider- 
ations applied to Figure 3.2b we obtain 


and from (3.3) we obtain 


C Q sinttp = h 
C sinQ!= h 


(3.5) 


Therefore 



(3.6) 


(3.7) 


and for small angle of attack 


a = 



cosA 


(3.8) 


The induced drag measured in the flight direction and the corres- 
ponding one measured in the direction normal to the leading edge are 
not the same. 


The lift is the same in both cases, but the downwash angles are not. 

The downwash velocity at the lifting line [see 2.3.1 and 2.4.4] 
must be computed with respect to the flight direction, since it is the 
wing span measured with respect to such direction that will determine 
the downwash velocity. 

An easy mistake would be to extend the use of the well known 
expression for the downwash angle (for an elliptic lift distribution) 
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(3,9) 


W 

V r 


0 


0 

TT A5L 


0 0 
to the case normal to the leading edge 

W C L 


v tt m. 

2 

In fact, using equation 3.4 and since Al = fRg / cos A , we find 


that 


J 0 


rrfR r 


Tt8L r 


(3.10) 


0 "“ v 0 

This implies that the induced drag measured in the two directions 
is identical, since 


and 


d. - x-2- 

1 o n *o 


d i‘ 


= & 


rr5t 


where 


d. = section induced drag, flight direction 

1 o 

d^ = section induced drag, direction normal to the leading 
edge. 

The downwash velocity at the lifting line is induced by 
trailing vortices aligned with the flight direction. 

We shall call the downwash measured this way W^ . 


The downwash angles in the desired direction can bow be computed 

by dividing by the corresponding velocity. Thus, the section 

induced drag measured in the flight direction is given by 

W. 


- £ 


V, 
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and, in a direction normal to the leading edge 


or 



W. 
x 

V q cosA 


d. 


(3.11) 


d. = V 

x cosA 


(3.12) 


Equation 3.12 shows that the induced drag measured normally to the 
leading edge is greater than the component in the flight direction. 


In a symmetric swept back (or swept forward) wing, the two side- 
components would cancel each other so that the total induced drag actu- 
ally experienced by the wing coincides with the component in the flight 
direction. For the oblique wing, the lack of symmetry introduces a side 
component of the induced drag whose section magnitude is given by 


(d.) = d. sinA 

v l y i 


d . tanA 
0 


(3.13) 


Therefore, for the case of the left wing forward, the side force 
contribution due to the induced drag of the wing is 


(F ) = “ ( D .) 0 tanA ( 3 . 14 ) 

In addition to the side force (F^^ due to the drag component in 
the y direction, a second term must be introduced when the wing is 
swept and at an angle of attack. The wing rotates perpendicularly to 
x^ and, therefore, parallel to (Fig. 3.3). 

Let us now consider the wing position w.r.t. the stability axes for 

the case a 0 and A i- 0 (Fig. 3.3). This position can be better 

visualized by considering these two steps: 1) a rotation of the aircraft 

about its y = y, axis by an angle a ; 2) a rotation of the wing 
so u 
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about z by an angle A . The final position of the wing is no longer 
parallel to the y g axis, the left (forward) tip being higher than the 
aft one. By idealizing the wing with a straight line and considering 
its projection y^ onto the z g , y g plane, we define the angle y^ 
(Fig. 3.3b), The lift vector being by definition perpendicular to the 
wing axis and to the velocity (parallel to x g ) , will therefore be 
banked by an angle producing a contribution (positive in this case) to 
the side force. This feature would not be present in a flying wing 
since, by banking the wing by y^ , it would be possible to realign lift 
and gravity, whereas in a conventional configuration the lift produced 
by the tail and the different inertia properties produce a somewhat more 
complex picture. 

We can now quantify this side force for the case of a straight 
rigid wing, by introducing some geometric considerations. 

Let us consider the wing for Q!^ = A = 0 and assume for simplicity 

that our stability axes have their origin at the wing pivot; in this 

case the y axes would be aligned with the wing span, 
s 

Let us now rotate a„ about y ; we obtain the set of axes x, , 

0 J s b ’ 

y^ , z^ , where y^ = y g is still aligned with the wing span. 

The rotation matrix for such transformation is 


r x " 

x b 


cosQIg 

0 

- sinQ 


X 

s 

y b 

= 

0 

1 

0 


y s 

z b 


sinQL 
L 0 

0 

COSQ^ 


z 

s _ 
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a) 


i 


Wing position for A = 0 

= 1 ===== wing position for a = 0 

and A^O (left wing forward) 



Figure 3.3 - Side Force due to Wing Rotation 



Figure 3.4 - Effects of Wing Bending on Side Force 



w 


' w 


If we now rotate the wing by A about its pivot, the wing axes 
will be defined by the following matrix transformation 

cos A sinA 0 
- sinA cosA 0 

0 0 1 | 1 z, | | sina„ 0 


w J l 



*b 



y b 

= 


_ Z b . 



0 


cosA 

sinA 

- sina Q cosA 

cosO!^ 

cosA 

sinA sinQ^ 


0 


cos a, 


'0 


2 

L s 


X 


X 

w 


s 

y w 

= [R] 

y s 

z 


z 

w _ 


L s _ 


The lift vector is perpendicular to the plane defined by and 

y^ and in terms of the stability axes 


U 0 = U 0 


(^- unit vector) 


■y^ (in stability axes) = [R] 


0 

1 

,0J 


m 


sinA cosQ 


co 


sA 


'0 


sinA sina. 


O' 


Let us define 


L = 


L 'zL 


where 'z'^ is a. unit vector 
3 


to 6? and ‘y’ 
s 


or 


•z 1 . • x =0 
3 s 


/S n 

z • y = 0 
3 w 
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since “S' (in stability axes) = 


I" 


0 j and 'z > 3 = 

\o! 


x 

z^ I we have 

y 


V 


z 3 ' ° 

X 


cosA z 2 + z 3 sinA sina^ = 0 

y z 


therefore 




(0 

- tanA sin a, 

u 


01 JT TT . 2^. 
\ 1 + tan A sin 


We can now define 'y' to "x 1 and 


% ■ \ 


'3=0 

x 


*3 • *3 - ° 


- y 3 tanA sina fl + y ^ 

y z 


= 0 


/ 0 


ft- 1 


V , 

1 + tan A sin QL 

\iau/isfiuu.( 0 


£„ = £ =10 
3 s 

\ 0 1 


The aerodynamic force in the x^ , y^ , z^ reference axes has the 
following components 
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We can now express these in terms of our stability axes 



0 

1 

2 2 

1 + tan A sin 

tanA sinCi^ 

2 2 

1+ tan A sin ct 
0 


0 

- tanA sinGl^ 

2 2 

1 + tan A sin a 

1 

2 2 

1 + tan A s in 



The value of y w.r.t. the stability axes x g 

to be 



z is found 
s 


tanA s inCi. 

sin y = j ? (3-15) 

(1 4- tan A sin a^) 

and for small angles 

y Q = a Q tanA (3.16) 

The corresponding lift contribution to the side force is 

<y L = (L) Tot Sin ^0 ~ (L) Tot a 0 tanA (3 * 17) 

All the previous analyses can be extended to the case of a. flexible 
wing having a built-in dihedral (and neglecting the twist contribution) 
by superposing a correction to the angle (Fig. 3.4) 

y(y) = y q + <p(y) cosA (+ bending modes) (3.18) 
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where 


<p(y) “ geometric dihedral at A * 0 

In this case the side force has to be evaluated by means of an 

integration along y g . The equilibrium condition may be reached by 

banking the aircraft by an angle <£q which, in terms of stability axes, 

corresponds to introducing a gravity component in the y direction. 

s 

For a straight rigid wing with no twist, considering total side force 
(contribution of the oblique wing itself and contribution of the geome- 
try of the wing rotation) and gravity, the equilibrium equation is 


Y = (F y ) D + (F y ) L + mg sin * <L) Tot « 0 tanA + (F y > D + mg sin <& Q 


4 >, 


0 mg 


<«Tot “0 tanA+ <Vd 


Since for the trimmed condition (L)_ = mg 


6 = 


a Q tanA + (F y ) D 


3.1.3 Wing Rotation. 

Figure 3.5 shows the geometry of the wing rotation about the pivot 
and the moment arms of the aerodynamic forces with respect to the center 
of mass (C.G.) of the aircraft. 

The aerodynamic force, as assumed in this study, is applied at the 
quarter chord. Since the wing is straight we can also assume that the 
line joining all quarter-chord points is straight. Such line is repre- 
sented by line "m" in the unswept condition, and by line "n" in the 
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swept one. The moment arm x about x^ la used when computing yawing 
and rolling moments as well as when computing the section perturbation 
velocity due to rate of pitch and/or rate of yaw. 

The moment arm y , about the y g axis, is used when computing 
angular velocities due to rate of pitch and the pitching moments. 

In determining the pitch stiffness C [Ref. 21] it is convenient 

"a 

to refer to the quarter-chord line for the unswept case, the moment arm 
about such line is indicated by y^, . 

Since the wing rotation affects the magnitude of the above moment 
arms, it will also modify the magnitude of all aerodynamic moments. 

The wing rotation also affects the position of the ailerons with 
respect to the aircraft centerline. Figure 3.6 shows the ailerons and 
section of wing which is affected by their deflection. Figure 3.6a 
shows the case for zero sweep. The shadowed area is defined by two 
straight lines aligned with the free stream velocity V,^ passing 
through the aileron's outer and inner stations. When the wing is swept, 
(Fig. 3.6b) we can notice that the two straight lines defining the area 
affected by the aileron are now passing to the left of the corresponding 
point on the quarter chord line for the unswept case. 

The result of this is an increase in the moment arm about the 
longitudinal axis for the left aileron, and a decrease in the moment arm 
for the right wing. The position of the pivot behind the quarter chord 
line will reduce and eventually eliminate such an effect. 

By simple geometric considerations, we can compute the quantities 
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X 



X 

s 


y s = 


x 


0 


y 0 = 


c = 

P 

d 

X = 

y 

y MAC 

Geometric 








stability axes 

axes parallel to stability axes centered at wing root 
quarter- chord, 

pivot-C.G. distance 

pivot quarter- chord distance 

moment arm about x 

s 

moment arm about y 

s 

moment arm about m (quarter -chord line for A = 0) 
expressions for the moment arms: 
x - y Q tanA + d sinA 
y = c + d cosA - y n tanA 

J p 0 

y MAC = V d(c0sA " 1} " y 0 tanA 


Figure 3.5 - Geometry of Wing Rotation and Moment Arms 


51 





52 




Figure 3.6 Influence of Wing Rotation on Ailerons Geometry 



and 


and C 2 , and (C^ and being the same as 

for the case when the ailerons are symmetric with respect to the 
wing centerline). 

3.2 Derivatives 

3.2.1 General Methodology. 

The perturbation quantities considered in a stability analysis are: 

- velocity in the x direction: u 

- velocity (sideslip) in the y direction: v 

- velocity in the z direction: w 

. ^ 

- rate of roll (or angular velocity about the x axis): p 

- rate of pitch (or angular velocity about the y axis): q 

- rate of yaw (or angular velocity about the z axis): r . 


Of these quantities, the second and third are generally normalized 
with respect to the free stream velocity . In doing so the new 

perturbation quantities are 


a t_^_ p ^ 

V H V 

TO TO 


(3.20) 


which correspond respectively to a change in the angle of attack and a 
change in the sideslip angle as can be seen in Figure 3.7. 


While the effects of a perturbation u , (3 , or r must be consid- 
ered separately, the perturbations p and q can be dealt with in the 
same fashion as for a , since in fact they all affect the local angle of 
attack, as will be shown in the next section. 


We shall now describe the general methodology for evaluating the 
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stability derivatives in a systematic way according to the previous 
subdivision; a much more complete discussion can be found in References 
15, 21, and 22. Such methodology applies to dimensional stability 
derivatives. 

We refer to Appendix A for the relationship between non-dimensional 
and dimensional derivatives, and to Appendix B far detailed calculations 
of the non-dimensional stability derivatives. 

3.2.2 a , p , q Derivatives. 

Let us now concentrate on a section of our wing. Whenever we 
consider a positive perturbation CC in the local angle of attack, the 
wing actually experiences a new perturbation velocity w . The result 
of this velocity w is that the free stream velocity becomes V having 
a new direction which differs by an angle a from the previous one. 
Consequently, the lift and drag components of the aerodynamic force must 
now be referred to this new direction. This result is shown in Figure 



Figure 3.8 Section lift and drag force for an a perturbation. 
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I 


Let us now evaluate the forces with respect to the stability axes 
system for this perturbed case. 


X ~ i sina - d cosa ~ i a - d 
Y = Ji y - d tanA 
Z -- - i cosa - d since ~ [1 + d a] 
and the corresponding moments are: 


(3.21a) 


Z x 


M s = " Z y MAC 

= -X x + Yy 


(3.21b) 


For the equilibrium case, the equations are 


Y o = Vo - d o tan(A) 

Z 0 ^0 (3.22) 

£ S = z o* 

M 0 S Z 0 y MAC 

N o “ x o x + Y o y 

The changes due to the perturbation a are therefore given by 


AX = X - X Q = i a - Ad 
AY - SLy - A y Q - Ad tanA 
Az = Li + da 
A^g = AZX 

- - ^ y MAC 


AN S = - AXx + AYy 


(3.23) 
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Because of the small perturbation assumption these equations can be 
linearized by expanding them in a Taylor aeries about the equilibrium 
condition (or CX= 0) and retaining only the first-order terms. For the 
first equation, this approach leads to 


CC (3.24) 

a- 0 . 

and the same can be done for the other equations. 

This is the classical assumption of linear aerodynamic theory which 
is to accept the following approximation for stability derivatives: 


Ax. = 


5 

da 


ia - A d 


A J&\ 
b \db / b= ( 


(3.25) 


where 


A = any aerodynamic force (or moment) 
b = any perturbation quantity. 


Consequently, the stability derivatives at our section for a per- 
turbation a are given by 


v = /M \ . Y ) . z =(M) 

a \3a/a=0 ’ a \ da Ja= o 1 a \3a/a= 


0 


(3.26) 


£ s = z a x 
b a “ 


M s z a y MAC ; N s_, x a X+Y a y 

\X 


where 


^aa-Ad) 


a= o 


l + M a- — 
Z + 3a da 


a- 0 " ^0 


\3a / a= 


(3.27) 
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Y a = 


55 <iy ' Vo ' ^ tanA) 


a= o 


(3.28) 


H y + * • H timA .a- 0 * l Jo + *o to ‘ Is tanA 


z a = 


30 : 


(i + dcc) 


a= 0 


Is + H a + d 


a= 0 - ^a + d o ( 3 * 29) 


3a 3a tanA+ c ® nst terms ^a= 0 = tan ^ 


(3.30) 


The same approach can be used when studying and perturba- 

tions. In both cases the perturbation will produce a change in the 
angle of attack. Such change will no longer be constant along the span, 
but will be 


for the pitch, and 


a = 


_SL_ 


v. 


TO 


(3.31) 



(3.32) 


for roll. 


Although the distribution in the angle of attack is antisymmetric, the 
corresponding increment in lift, and consequently drag, because of the 
shape of the FWID will not be asymmetric. The asymmetry in lift, case of 
a symmetric wing in roll, does not produce any variation to the lift 
vector. Loss of asymmetry introduces a change in the total vertical 
force Z as well as in X . 

The same approach used for the a derivatives can be applied to the 
evaluation of roll and pitch derivatives. 
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The only differences are: 


a) the section perturbation angle of attack which is no longer a 
constant as in the previous case, but varies according to Eqns. 
(3.31) and (3.32); 

3 

b) the linear operator must be replaced, for the pitch by 


_3_ _ _3_ 3a _ y 3 

3q 3a 3q " V T0 3a 


and for the roll by 


(3.33) 


3 _ _3_ 3a _ x 3 

3p 3a 3p V tQ 3a 

The angle y remains unchanged. 


(3.34) 


The wing aerodynamic derivatives with respect to the stability axes 
can be obtained by integrating the section values over the span. 


3.2.3 p Derivatives. 

The effect of introducing a sideslip angle (3 corresponds to a 
negative A A for the wing (Fig. 3.9). Consequently, for the wing 


_3_ 3_ 

3p 3 A 


(3.35) 


The evaluation of the p derivatives is not trivial when consider- 
ing the aerodynamic moments. In this case, because both FWLD as well as 
the distribution of lift due to twist will be affected, it is not as 
simple to derive analytic expressions as for the previous cases. A 
numerical calculation would be possible by considering results at two 
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different wing sweeps and then comparing them. For the cases where no 
moments are involved, the total aerodynamic forces considered are 


F = - D cosp - (f ) r + (f 1 
x \ y/L A y/w 


sing ^ - D - 


( F y)L + ( F y) 


W 


F = - D sing + 

y 


( F y)l + ( F y)» ] ^ + [( F y)L + ( F y)„ 


P 

(3.36) 


F Z = - < L) Tot 


In addition, the following approximation can be made. 


(L) 


Tot 


a) 


Tot 


A= o cos A 


(3.37) 


fD) = D + (L) Tot ^ C L ) Tot _ 

^ Tot ' D 0 + TT& °0 + 


^Tot ^L^Tot 


4 , 
cos A 

TT7R 


(3.38) 

2 

and, since 7R = (7R )^_ Q cos A , it is possible to evaluate the g 
derivatives of these forces. 


3.2.4 r Derivatives. 

The new yawing derivatives introduced by the oblique wing are due 
to two factors: a) side force and b) aerodynamic coupling between rolling 
and pitching moment. A rate of yaw increases the local stream velocity 
according to the relation 

AU= - x * r (3.39) 

Therefore, for the dynamic pressure the variation is 
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H.O.T. 


q(y) - 2 plj2 = 2 p(u o + Au)2 = 2 P K + 2u o Au+ 


q(y) - \ p(u 2 0 ■+ 2U 0 Au) 


(3.40) 

( 


Substituting (3.39) into (3.40), we obtain 


q(y) = q(l - 2 ^ r) (3.41) 

The linear antisymmetric change in q(y) due to the perturbation 
r will now affect the aerodynamic forces due to both the flat wing and 
twist (including effects of dihedral and camber) contributions. Since 
the antisymmetric variation in q(y) is applied to a symmetric (cruise 
condition) distribution, the variation in the loading distribution will 
also be antisymmetric, like for the case of a symmetric aircraft. There- 
fore, when first order variations only are considered, no changes will 
occur in the wing total aerodynamic forces, but all three aerodynamic 
moments will be affected. 


Since q(y) is the only term which changes in the expressions for 
X , Y , and Z , the evaluation of the yaw derivatives is reduced 
to computing 


ax = ^0 

dr q 3r 


x 


q|i ' 2 ^ r 


■ - x o • 2 u 


X 


0 


(3.42) 


and similarly for y , and Z , 


= - Y 2 ~ 
37 *o u 


(3.43) 


0 
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X 


( 3 . 44 ) 


c>Z 

S' ' - z o 2 u 


0 


A simple integration over the span of the last three expressions 
will determine the yaw derivatives of the aerodynamic forces. For the 
£ , M and N c moments, the expressions to be integrated are 

D J D 


f s = - z o 2 ^* 


M s ' + z o 2 T 0 y 


( 3 . 45 ) 


N S + x o 2 u Q x " y o 2 u Q y 
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3.3 Summary . 


In this Chapter we have introduced and evaluated the two side 
force contributions, one due to induced drag and one due to wing 
rotation, that an oblique wing aircraft would experience. The effects 
of wing rotation were then analyzed. The general methodology for 
computing stability derivatives was also discussed. 
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IV EQUATIONS OF MOTION 


4.1 Rigid Body 


The equations of motion for - the aircraft can be derived from Newton's 
Second Law of Motion, which states that the summation of all external 
forces acting on a body must be equal to the time rate of change of the 
momentum of the body, and the summation of the external moments acting 
on a body must be equal to the time rate of change of the moment of 
momentum (angular momentum). The time rates of change are all taken 
with respect to inertial space. These laws can be expressed by two 
vector equations. 


v-!£? 

m 

(4.1) 

I 

H =X)T 

(4.2) 


where I indicates the time rate of change with respect to inertial 
space. 

Now, the external forces and moments consist of equilibrium values 
plus a perturbation value which stems from a difference from this equi- 
librium condition. Thus, 


F = F q + AF 

t = t q + At 


(4.3) 


In the dynamic analyses to follow, the aircraft is always considered 
to be in equilibrium before a disturbance is introduced. Thus, Fq and 
Tq are identically zero. 
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The equilibrium forces consist of lift, drag, thrust, and gravity, 
and the equilibrium moments consist of moments resulting from the lift 
and drag generated by the various portions of the aircraft and the thrust. 


The following observations are made to achieve a linear analysis. 
First, the aircraft is in an equilibrium condition before perturbation. 
Second, the mass of the aircraft remains constant during any particular 
dynamic analysis. Third, it is assumed that the aircraft is a rigid 
body. Fourth, it is assumed that the earth is an inertial reference 
and unless otherwise stated, the atmosphere is assumed to be fixed with 
respect to the earth. The time rate of change of the velocity vector 
with respect to the earth, in the assumed reference axes, is given by 

I S 

V = V + tt) X V (4.4) 


where S indicates time rate of change with respect to the assumed 
reference axes system (non inertial); to is the angular velocity of the 
aircraft with respect to the earth; and x signifies the cross product. 

Similarly ^ s 

H = H + to x H (4.5) 


V and w can be written as 

w = w. + "Aw 

v = v Q + Av 


(4.6) 


where Wq and V^ represent the steady state conditions and Aw and 
AV are the perturbation quantities. 


Their components in stability axes are: 
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The angular momentum is given by the dot product of the inertia 
tensor, I , and the angular velocity vector, U) 


H - I • u> 


(4.8) 


and Equation (2.5), 1 being time invariant because of assumption 

number four, can be rewritten as 


I 

H = I 


S .•¥ 

— + CD X (I 


OJ) 


(4.9) 


The inertia matrix of an oblique wing aircraft differs from the one 
of a symmetric aircraft because of the non-zero product of inertia I . 
For an oblique wing aircraft, the inertia matrix in body axes is 



Because our stability analysis is based on the stability axes system 
previously introduced, it is necessary to transform the inertia matrix 
from body to stability axes. This transformation is carried out in detail 
in Appendix E. 

Because of the transformation the inertia matrix in stability axes 
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will no longer have a zero I . In fact, as shown in Appendix C, it 

yz 

is given by 


I' = - I sina 
yz xy 


(4.10) 


where the prime denotes stability axes. 


The equilibrium condition we are going to consider in our analysis 
is that of fl4traight level flight, therefore 


P = Q = R=V 0 = W 0 =0 


(4.11) 


For this condition, after substituting (2.4) into (2.1), and after 
some algebra, we obtain 


if-i F 

m my 


v + u„ r + 


w - q u Q> 
linear terms 


qw - rv 


ru - pw 


vp - qu 


non linear terms 


(4.12) 


Similarly, Equation (2.2) becomes 


I' $ + 1' q+I' r 
xx xy xz 


I ' p + I ' q+I' r 
xz yz zz 


linear terms 


(I* q - I' r) q + (I' - I' )qr 

x z xy zz yy 


M = I* p + I’ q + I’ r + I' (r 2 -p 2 )+I’ qr+ (I ' -I’ )pr 

xy yy yz x z xy xx z z 


I' (p 2 - q 2 ) - I' qr+ (I' -I' )pq 

xy sr n xz ^ yy X.X. 


non linear terms 


(4.13) 


It is now necessary to expand the applied forces and moments and to 
express them in terms of the changes in the forces and moments that cause 
or result from these perturbations. These latter forces are usually of 
an aerodynamic and gravitational origin. 


If it is assumed that, as the disturbances are small, the partial 


69 



derivatives are linear, the differentials can be replaced by the actual 
increments , and F and T can be written as 


F - AF =• 


M- ) h 
x \ 3b / i 

V=o 

i 


(4.14) 


t = At * 2). 


3t 


where 


3f 

lSb i/b.= 0 
x 


and 


j 3t \ 

\ ab i/b = 


i \ 3b. I i 

X b ,=0 
x 


represent the aerodynamic stability 


derivatives evaluated at the steady state condition when the perturbation 


variable b^ is zero. 


The linearized set of equations therefore is; 


u= E 


/ 3F 


x 


x 3b./ i 

x b .=0 
1 


*+ u 0 r 


1 'b.=0 

X 


(4.15) 


3f 


W - qu. 


- ^ s? K 


b.=0 

x 


I' p + I' q + I' r = 
xx xy xz 


Si 


^ -i 

' x /b.=0 

x 


70 



I' t> + I' q + I' r = £. /— )b 

vtt 4 trir * ir*7 i 


xy yy yz 


i \3b. /"i 
ly h .=0 

i 


-A 


I' p + I' q + I ' r= 
XZ yz zz 


Si £ - 


/ 1 

ly b.=0 

i 


(4.15) 


In the set of Equations (4.15) the terms q and I^_ p } not 

present in the case of a symmetric aircraft, introduce the inertia 

coupling between pitch and roll motions. The terms I 1 q and I' r 

yz y z 

have also a coupling effect and seem to make the coupling between 
longitudinal and lateral dynamics even stronger, but their coupling is 
only apparent since it depends on our choice of axes system. 


In addition to the inertia coupling existing in symmetric wing air- 
craft between roll and yaw, oblique wing aircraft experience an inertia 
coupling pitch and yaw motions. In table 4.1 we report the matrix 
representation of the set of Equations (4,15) in terms of Laplace trans- 
form. 


4.2 Effects of Flexibility. 

The analysis done so far applies to a rigid aircraft. 

Now, it is known that the stability and control characteristics of 
flight vehicles may be profoundly influenced by the elastic distortions 
of the structure under aerodynamic load. Many of the important effects 
of distortion can be accounted for simply by altering the aerodynamic 
derivatives. The assumption is made that the changes in aerodynamic 
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TABLE 4,1 Laplace Transform of the Equation of Motion. 


loading take place so slowly that the structure is at all times in static 
equilibrium. This is equivalent to assuming that the natural frequencies 
of vibration of the structure are much higher than the natural frequen- 
cies of the rigid-body motions. Thus, a change in load produces a 
proportional change in the shape of the vehicle, which in turn influences 
the load. 

When the separation in frequency between the elastic degrees of 
freedom and the rigid-body motions is not large, then significant inertial 
coupling can occur between the two. In that case, a dynamic analysis is 
required, which takes account of the time dependence of the elastic 
motions. 

The elastic motions have no inertial coupling with the rigid-body 
motions except through I . However, it has already been assumed, in the 
previous paragraph, that such time rate is second-order and negligible 
in the small-perturbation theory. 

The only coupling existing between elastic and rigid body motions is 
due to the aerodynamics. 

The aerodynamic derivatives associated with the deformations of the 
airplane are of two kinds: those that appear in the rigid body equations, 
and those that appear in the added equations of the elastic degrees of 
freedom. 

In our analysis we shall use a quasi-steady approach neglecting 
unsteady aerodynamic effects and assuming that only the wing is flexible. 

Therefore, we shall also neglect the added equations of the elastic 
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degrees of freedom and refer to Ashley, E'tkin, and Blakelock [15, 21, 22] 
for the complete study. 

In a quasi-steady analysis the effects of flexibility on the wing 
are basically reduced to static bending and torsion. 

Because of this bending, the aerodynamic center shifts forward 
changing the value of all the aerodynamic moments and their related 
derivatives. The side force due to lift is also affected. These static 
effects have been included in our numerical analysis and will be 
described in that section. 
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V NUMERICAL ANALYSIS 



In the previous chapters we have discussed qualitatively the influ- 
ence of an oblique wing on the aerodynamics and dynamics of aircraft. 

In this chapter we shall attempt to quantify these new properties and 
evaluate their influence on the stability of an oblique wing aircraft. 


The aircraft configuration that has been studied is the Boeing 
single-body yawed wing aircraft model 5.3 (Fig. 5.1), that was studied 
for NASA [Ref, 24], Unlike what is shown in Figure 5.1, our analysis, 
as well as the one done by Boeing, was carried out for the 'left wing 
forward" case. The nominal configuration and flight condition for the 
simulation were: 


Mach number =0.8 

Altitude = 20,000 ft (6096 m) 

Gross weight = 400,000 lb (181,440 Kg.) 

Wing sweep = 45 degrees 

C.G. location = .355 (M.A.C.). A (body station 57.8 m) 

TrO 

The inertia properties being referred to body axes were transformed 
to the equivalent in stability axes according to the transformation 
described in Appendix C. 
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The .355 M.A.C. center of gravity location in the nominal con- 
figuration produced a longitudinally unstable vehicle. This con- 
dition was selected by Boeing in order to achieve a more satisfactory 
overall vehicle design; since the design included a longitudinal 
S.A.S. this was also the case in the Boeing SST design in the subsonic 
regime. 

In the unstable (without SAS) configuration, the normally 
complex short period roots migrate significantly; one moving to the 
positive real axis and one combining with the "rolling convergence" 
root to form a new complex pair in the LHP. 

In order to study the effect of the oblique wing under more 
"normal" conditions, we chose to modify the nominal configuration 
so as to achieve a stable configuration for the zero sweep ( A = 0°) 
condition. This was accomplished by moving the C.G. forward. 
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5.2 Aerodynamic Results 


In chapter II we described two different methods for evaluating the 
spanwise lift distribution for the oblique wing and how to obtain the 
spanwise induced drag distribution from the knowledge of the lift. 

Two numerical programs, based on these two methods, have been 
written: OWSD/ST computes the stability derivatives for a rigid oblique 
wing aircraft by means of strip theory; OWSD/LT , instead, applies the 
linear theory approach and can also include the static effects of a 
flexible wing. 

We shall now report the numerical results obtained with the two 
methods, for a wing at a skew angle A of 45°. 

In our strip theory analysis, the wing was assumed to be rigid; 
therefore, the results obtained with that analysis will be compared 
against the rigid wing case of linear theory. The, we shall compare 
the results given by linear theory for the rigid and flexible wing. 

5.2.1 Comparison between Strip Theory and Linear Theory (Rigid Case). 

Figure 5.2 shows the spanwise lift distribution at A = 45° as 
computed by means of the lifting line method, for the flat wing as 
well as for the case including built-in twist. The twist is linear, 
and the tip values were chosen so that the resultant lift was acting 
at the wing centerline. In a real design the airfoil camber distribution 
would be so that the cruise lift is acting at 50$ chord. The assumed 
twist distribution certainly is not the optimal one, leading to a 
minimum induced drag, but, at least, it produces no rolling moment 
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and, for the purposes of our analysis, is a satisfactory cruise 
condition. 

The twist values, for the unswept wing, are 

Left tip = + 2.6 degrees 

Right tip = - 3.35 degrees 

The downwash angles for the flat wing computed by the two methods 
are shown in Figure 5.3. The lifting line approach (Fig. 5.3a) shows 
a large variation in the spanwise distribution of the downwash; it also 
indicates that the right wing, beyond the 60 j half span, experiences an 
upwash which, in terms of horizontal force, results in a thrust. 

The induced drag distribution is shown in Figure 5.5 where the 
negative value corresponds to the thrust mentioned before. 

The results obtained by means of strip theory show a slight varia- 
tion in the downwash angle distribution and the corresponding induced 
drag distribution appears to be almost symmetric about the wing center- 
line. The drag distribution, as obtained with the lifting line theory, 
produces a yawing moment which tends to unskew the wing. This effect 
is not shown in the results obtained with strip theory. 

The tendency to unskew the wing has been experienced also during 
wind tunnel tests run at NASA-Ames Research Center. Since the model 
was not a pressure one, it is not possible to tell whether this was 
due only to a higher drag on the forward wing, or also to a thrust 
force in the aft one. 

The side force distribution (Fig. 5.4) predicted by the two methods 
differs since the drag terms differ. The side force predicted by the 
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lifting line theory increases the unyawing moment produced by the drag, 
instead, the strip theory result shows the opposite effect. 

We can therefore conclude that the spanwise induced drag distribu- 
tion, as computed by means of strip theory, is a poor approximation. 

We shall now continue in the comparison between the two methods 
for the cases when perturbations are considered. Unfortunately, the 
numerical program computing stability derivatives according to strip 
theory computes only the nondimens ional lift distribution. Thus it 
is only possible to compare the results obtained with the two methods 
qualitatively. 

The first perturbation quantity considered is Ct . 

Figures 5.6 and 5.7 show respectively the A F x (positive forward) 
and AF • For the sake of clarity we recall that, for the perturbed 

y 

case AF x is given by 

Af^ ^ L q a - D 

Figure 5.8 shows the change in lift due to a roll perturbation. 

According to strip theory, no change occurs in the lift at the air- 
craft centerline, whereas the lifting line solution shows, for positive 
rate of roll, a negative AZ at the aircraft centerline. 

A greater difference shows up in the change in lift due to a pitch 
perturba.tion. This greater difference is also due to the fact that, in 
the case of strip theory, the vertical velocity introducing the aero- 
dynamic twist is computed at the quarter-chord, whereas, in the lifting 
line method, it is evaluated at the 3/4- chord point. 
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Figure 5.10 shows the total lift distribution for a negative side- 
slip of 5 degrees. 

We shall return on the sideslip case in the next section, where 
we discuss in detail the behavior of the oblique wing in sideslip. For 
the moment, we shall only point out how good the agreement is between the 
lifting line theory results and those obtained with the empirical method 
described in section 2.5. 


5.2.2 Comparison between Rigid and Elastic Wing (Linear Theory). 

The lifting line method described in section 2.4 implies the know- 
ledge of the station 2-D lift slope. For this reason a value for this is 
considered an input data for the numerical program OWSD/LT. In our analysis 
the section liftslope has been assumed to be constant along the span. 

This is not true in most cases and specially for our model, since the 
thickness to chord ratio varies along the span. The value of the 2-D 
lift slope has been adjusted so that the wing lift slope coincides with 
the value computed by Boeing [Ref. 24] in its study. No built-in 
dihedral was assumed, but only a linear twist was introduced, satisfying 
(as mentioned in the previous section) the zero rolling moment condition. 

In our model, the pivot location (50$ root chord) does not coincide 
with the quarter-chord line (where the lift force is assumed to be acting); 
therefore, when the wing is skewed, the wing center line does not coincide 
with the aircraft centerline. For that reason the right wing has a span 
larger than the left one. Therefore, the lift distribution for the cruise 


88 



condition, satisfying the equilibrium condition about the roll axis, can 
no longer be symmetric in order to balance out this asymmetry in the two 
semiwings. This asymmetry can be noticed in Figure 5.11 which shows the 
spanwise lift distribution at A = 45°(rigid and flexible wing)for this twist 
distribution and for an angle of attack of 3.75 degrees (.0652 rad.). 

Since the twist was computed for the rigid wing case, the lift distri- 
bution of the elastic wing will produce a strong positive rolling moment. 

The flat wing case is shown in the next figure (Fig. 5.12) for the same 
angle of attack. 

It is interesting to notice that the flexible wing case produces 
almost no rolling moment about the aircraft centerline (The resultant 
lift vector is applied at the quarter-chord station having a distance 
to the right of the longitudinal axis equal to 0.72 ft). 

Such a small difference does not even require a built-in twist, 
but only a little aileron correction. At this point it should be pointed 
out that the built-in twist is, for an oblique wing aircraft, a poor 
design solution. In fact, since the twist measured in the flight di- 
rection varies with the cosine of the sweep angle, the twist correction 
is maximal when the wing is unskewed ( and then we do not need any twist 
correction since the FWLD is symmetric) and decreases when we skew the 
wing. 

A built-in dihedral is, therefore, a far better solution, producing 
no twist, for the unskewed case, and increasing the twist correction 
when the skew angle increases. 

We shall complete the analysis of the cruise condition, as from our 
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numerical model, with Figure 5.13. Figure 5.13a shows the downwash 
distribution for the flat wing, where the negative values correspond 
to the upwash. The flat wing drag distribution corresponding to this 
downwash is shown in Figure 5.13b. The flexible wing shows a much higher 
drag on the forward (left) wing and a thrust over a larger portion of the 
aft wing. This is expected; in fact the forward wing, by bending, in- 
creases the angle of attack and consequently the lift (this behavior 
is well known as the divergence problem of a swept forward wing), where- 
as the opposite happens on the aft one. 

The direct consequence of the higher lift of the forward wing is an 
increase in the upwash velocities induced on the aft one. Instead, for 
the aft wing, the decrease in lift will result in a decrease in the up- 
wash field induced on the forward one and, consequently, lead toward 
a further increase in the downwash of the aft wing. These are the rea- 
sons why the flexible wing experiences a higher downwash value on the 
forward wing, and a larger portion of the right one is affected by the 
upwash velocities. 

Then, when we multiply the lift (a positive value) times the corre- 
sponding downwash, in order to obtain the spanwise induced drag distri- 
bution, we multiply the already greater downwash times a larger value 
of lift for the left wing. For the right wing, instead, the lift magni- 
tude is smaller and this is the reason why, though affecting a larger 
portion of the aft wing, the thrust reaches a lower maximum value than 
the corresponding one for the rigid case. 

We shall now continue our comparison analyzing the results for the 
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cases simulating perturbed states as are assumed in our computation of 
the stability derivatives. 

The numerical program OWSD/LT evaluates the stability derivatives 
by first computing the distribution of the aerodynamic forces (lift, 
drag, and side force) for the perturbed condition, then evaluating the 
increment from the corresponding cruise distribution, and finally inte- 
grating these increments as well as their product with the corresponding 
moment arms in order to obtain the incremental moments over the entire 
wing span. 

The desired values of the stability derivatives are then obtained 
by dividing the results of the integrations by the perturbation quantity 
and by the appropriate nondimens ionalizing factors [See Appendices A and 
B], Then, derivatives are computed with standard expressions [Appendix 
B and/or Ref. 15, 21] and the yaw derivatives can be evaluated according 
to strip theory as described in section 3.2.4. Therefore, for both of 
them it is necessary to simulate the perturbed condition. Consequently, 
the conditions to be simulated have been reduced to: angle of attack, 
roll, pitch, and sideslip. 

The CC - perturbation has been simulated by increasing the cruise 
angle of attack by 3 degrees (.05236 rad.). The new lift distribution, 
including twist contribution, is shown in Figure 5.14, 

The rigid wing shows a visible increase in the aft wing l±f$, due to the 
increase in the flat wing lift. For the elastic one, this effect is much 
less pronounced. These two behaviors can be predicted by looking at the 
shape of the FWID (Fig. 5.12); in fact, a change in a will introduce a 
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first order variation in the lift, proportional to these shapes. The 
variation, AF^ , in the horizontal force (Fig. 5.15a) is strongly af- 
fected by the lift contribution in that direction during the perturba- 
tion, as we have discussed in section 3.2.2; in fact 

Af = H a- (d - d rt ) 
x 0 

The increase in lift for the rigid wing is more pronounced in the 
aft part where it contributes to the thrust already experienced because 
of the upwash field, but on the forward wing its increase is not suffi- 
cient to offset the drag increase. 

For the elastic wing, the FWLD has a' more symmetric shape than for 
the rigid one; we would therefore expect the transition point from 
negative (drag) to positive (thrust) X-force to occur further outboard 
on the left wing, instead, the two shapes behave similarly, except for 
the magnitude. This is not the case, according to our results, and the 
reason is probably because of the downwash distribution: stronger down- 
wash on the forward wing, and consequently a higher induced drag to 
counteract the higher increase in lift; but also stronger upwash on the 
aft one with higher resultant thrust where the lift increment is lower 
than the rigid case. The same reasoning can be applied to the variation 
in the sideforce (Fig. 5.15b), where the two shapes look similar, al- 
though they differ in magnitude. 

The perturbations due to roll and pitch are simulated by intro- 
ducing an aerodynamic twist corresponding to a linear variation in the 
free stream vertical velocity [See section 3.2.2], The roll angular 
velocity is chosen to be such that the twist is equal to -5 degrees- at the 
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left tip and zero at the aircraft centerline [see Eq. 3.31], 

For the pitch velocity, the left tip condition remains - 5 degrees 
and zero twist now coincides with the wing station corresponding to the 
intersection of the quarter-chord line with the lateral axis y . 

The three different twists used in our study are shown in Figure 

5.16. 

Figure 5.16a shows the built-in twist; the aerodynamic twist corres- 
ponding to our roll simulation is shown in Figure 5.16b. The aerodynamic 
twist due to the pitch simulation (Fig. 5.16c) is not perfectly linear 
because the downwash velocities have been measured at the 3/4-chord 
points, which do not lie on a straight line. As we would expect, the 
aerodynamic twist introduced by a positive rolling motion increases the 
lift distribution on the aft wing and decreases it on the forward one. 


Figure 5.17 shows the spanwise lift distribution corresponding to 
the cruise angle of attack plus a twist given by the superposition of 
twists a) and b) of Figure 5.16. The increment in the vertical force is 
of greatest interest. Although such an increment has a contribution due 
to the drag, according to [see section 3.3.2] 


Az = + da 

where 


a = 


V 

TO 


such a contribution is negligible compared to the actual lift changes, 
therefore the AZ shown in Figure 5.18 can be thought of as the actual 
change in lift without any serious loss in accuracy. 
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From those distributions we notice that the rigid wing shows a 
positive variation on the aft wing, larger, in magnitude, than the 
negative one. The opposite is true for the flexible wing. Moreover, 
if we look at the shape of the corresponding case, as simulated by strip 
theory, we can observe how erroneous the strip theory approximation is 
for the flexible wing case. 

The variations Ax and Ay look quite complex (Fig. 5.19 a) and 

b)). 

The most important derivative computed from these two distributions 
is the yawing moment. In a symmetric aircraft the side force is always 
negligible. This, as we have proposed in section 3.1.2, is not true for 
an oblique wing aircraft, in fact the side force has a magnitude compa- 
rable with the horizontal one even in a perturbed condition. By inspec- 
tion, we can also notice that the yawing moment produced by the Ax 
and Ay distributions is, in this case, negative , and the rigid wing 
shows a larger magnitude. 

The same analysis can be extended to the pitch perturbation; there 
is only one difference: the way the aerodynamic twist is measured in our 
linear theory model. 

Eqn. 3.31 implies a linear variation in the aerodynamic theory; this 
is a common assumption and was used also in our strip theory analysis. 

The linear theory, as we mentioned earlier, implies the evaluation of the 
vertical velocities, induced by the angular velocity q , at 3/4-chord. 
Figures 5.20 and 5.21 respectively show the spanwise lift distribution 
corresponding to the cruise angle of attack plus the superposition of 
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the twists of Figures 5.16a and 5.16c. 


Figure 5.22 shows the A'Z distribution and 5.22b the AY distri- 
bution. 


In our model, the sideslip is simulated by increasing the skew angle 
by 5 degrees. This, as we have already mentioned, corresponds to a nega- 
tive sideslip according to the reference axes system chosen. Thus, when 
looking at the variations in the aerodynamic forces, we must change the 
sign when considering the sideslip case. Figure 5.23 shows the total 
lift distribution for the new skew angle. 

Increasing the skew angle, the angle of attack measured in the flight 
direction decreases (Eq. 3.8). In our case, the relationship between 
perturbed and cruise angle of attack is: 


a. 


50 ° 


- a 45 o 


cos 45 


cos 5 0 

The results shown in Figure 5.24 are of great interest. 


Increasing the skew angle (which corresponds to a negative slip), 
the total lift decreases and the spanwise distribution of such variation 
is asymmetric (Fig. 5.24). This asymmetry will result in negative 
pitching and rolling moments; therefore, a positive sideslip introduces 
positive rolling and pitching moments. This is an unstable behavior and 
can be visualized by considering the trend of the spanwise lift distribution 
when, the wing is skewed. 

The upwash field generated by the forward wing introduces the 
build-up in lift on the aft wing. 
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As we increase the skew angle, this build-up will shift more and 
more towards the aft wing tip. Thus, this shift will introduce negative 
pitching and rolling moments (left wing forward case). When we unskew 
the wing, the trend will be the opposite; consequently the lift build-up 
will shift towards the center decreasing the magnitude of the negative 
rolling and pitching moments. These are positive increments in the 
moments and the reduction in skew angle corresponds to the sideslip 
situation. The effects of flexibility worsen this undesired behavior. 

In our model, we have used linear built-in twist to stabilize the 
rolling moment. In this case also, we see how poor such a solution 
would be as compared to the case of a built-in dihedral. 

The effectiveness of the linear twist decreases with the cosine of 
the skew angle; thus, the amount of twist required at 45° to trim the 
aircraft would introduce a strong positive rolling moment when the wing 
is unskewed. 

Since a sideslip perturbation affects the total lift distribution, 
the contribution due to twist in a positive sideslip situation is an 
increase in the undesired positive rolling and pitching moments. Since 
the equivalent twist due to dihedral varies with the sine of the skew 
angle, the previous unfavorable situation now becomes a favorable one; 
in fact, a decrease in sweep will introduce a variation in the equivalent 
twist which results in a decrease in the forward wing lift and an in- 
crease in the aft wing lift, which is a stabilizing trend. The yawing 
moment due to the variation in the longitudinal force, AX. , is positive 
for our negative sideslip simulation (Fig. 5.25a), and this is a favor- 
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able result. 


From the shape of the side force (fig. 5.25b) we can observe that 
the contributions of flat and rigid wings to the yawing moment can differ 
even in sign. In fact, for the negative sideslip case, the rigid wing 
shows a decrease in side force whose resultant force is applied to some 
point on the aft (right) wing. The rigid wing, instead, shows the 
resultant side force acting on the forward (left) wing. Therefore, the 
rigid wing, in sideslip, experiences a side force that produces a desta- 
bilizing yawing moment (positive for the negative sideslip case shown in 
Figure 5.25b), whereas the elastic wing shows a tendency to produce a 
favorable jf&wihg’ moment . 

The effect of the wing rotation on the ailerons geometry has already 
been discussed in section 3.1.3. 

Table 5-Ia shows the ailerons dimensions for the unswept case, and 
Table 5-Ib shows the corresponding ones assumed in the computer simula- 
tion for the 45° skew angle; the quantities used in this Table are 
defined in Figure 3.6. 

Three aileron deflections are considered: 

1) left aileron deflected 5 degrees 

2) right aileron deflected 5 degrees 

3) both ailerons deflected; left 5 degrees, right -5 degrees. 

The elastic wing case considered so far assumed the elastic axis 
(E.A.) to coincide with the wing quarter-chord line. In addition to 
this elastic case and to the rigid one, we also considered, for the 
ailerons only, the case when the E.A. is at a distance from the quarter 
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YCNTOL = 93.63 ft CNTCOL = 2.2 ft 

Left Aileron 

YCNTIL = 73.93 ft CNTCIL = 3.66 ft 

YCNTOR = 93.63 ft CNTCOR =2.2 ft 

Right Aileron 

YCNTIR = 73.93 ft CNTCTR = 3.66 ft 


a) Ailerons Nominal Dimensions at A = 0 


Wing Station Distance from 

(non dimen.) A/C Centerline 

Left Aileron 
Right Aileron 

b) Computer Approximation at A = 45° 


TABLE 5-1 Ailerons Dimensions 
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chord. The latter represents the real situation, but, since at 45 skew 
angle the bending contribution to the wing twist is by far more important 
than the contribution due to torsion, the differences between the two 
elastic cases are small, at least as far as our analysis is concerned. 

We shall see how the same is true for the ailerons also. 

Figures 5.26 through 5.28 show the total lift distribution corre- 
sponding to the three aileron deflections considered, for both the rigid 
and elastic wing cases, having the E.A. coinciding with the quarter-chord 
line. 

The aileron effectiveness for these three cases is shown in Figures 
5.29a, b, and c. 

Once more, we can observe the peculiar behavior of the oblique wing. 
In fact, the forward aileron, affected by a stronger downwash, is less 
effective than the aft one, which is instead influenced by a strong 
upwash. The variations in horizontal force, Ax. , are shown in Figures 
5.30a, b, and c. 

Because of the lack of points in the areas of Interest, the shape of 
the curves looks unusual; nevertheless it is interesting to notice the 
AX variation introduced by the aft (right) aileron. 

The positive right aileron deflection produces an adverse yaw moment 
much smaller than the corresponding one produced by the left aileron 
(fig. 5.29 and 5.30). We have already noticed how the favorable upwash 
field decreases the drag on the aft wing; in case of right aileron de- 
flection the induced drag rise is not only lower than the corresponding 
one for the forward aileron, but the increment in lift, because of the 
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upwash, results in a thrust force outboard of the aft aileron. 

The same effect can be noticed in the case of antisymmetric deflec- 
tion (Fig. 5.30c). 

The negative deflection in the right (aft) aileron produces a 
decrease in drag, or equivalent thrust, followed by a drag increase on 
the outboard part of the wing. This behavior suggests the possibility 
of using the aft aileron only for control during the cruise condition; 
but this solution implies, first, the experimental confirmation of the 
analytic results and, second, that there are no adverse aeroelastic 
effects related to that solution. 

Figures 5.31 through 5.34 are referred to the elastic wing case 
having the E.A. passing through the wing pivot (50% root chord) and 
parallel to the quarter-chord line'. Very little difference can be noticed 
between these results and the corresponding ones for the other elastic 
wing case considered. 

In fact, the twist contribution due to bending is much larger than 
the contribution due to torque. 

5 . 3 Stability Derivatives and their Influence on the Natural Modes . 

At this point, the logic flow of our analysis would require the 
discussion of the numerical results of the stability derivatives. We 
shall postpone it to the next section and carry on, instead, a qualitative 
analysis of the influence of the stability derivatives on the natural 
mode of an oblique wing aircraft. 
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In Chapter III, the new stability derivatives due to the skewed wing 
were introduced without any attempt to define their importance to the 
dynamic stability and natural modes of the aircraft. This analysis had 
been carried out only on the basis of understanding the behavior of an 
oblique wing aircraft and its differences from a symmetric aircraft. We 
shall now discuss the result of a numerical investigation to determine 
the influence of the new stability derivatives, as well as of the usual 
ones, on the dynamic stability of an oblique wing aircraft. 

The set of six linear differential equations dervied in Chapter IV 
was numerically solved using an available computer program (GSA), which 
was originally developed by Lockheed Missiles and Space Company. The 
program solves a set of linear differential equations using the Laplace 
transform method and gives root locus, bode or time domain plots for the 
system. The characteristic polynomial is of the 8th order. 

The influence on the roots of the characteristic equation is now 
analyzed for each nonzero derivative by means of a root locus. In this 
way, it is possible to find the derivatives having influence on the 
natural modes. This study has been carried out for the new as well as 
the conventional derivatives, by varying one derivative at a. time, start- 
ing from zero to a value double that of the corresponding one reported 
in Table 5 -III, column 1. 

It is obvious that this is actually not possible in a realistic 
analysis, since the parameters affecting one derivative may affect other 
derivatives also; the range itself does not reflect a real case except 
for some of the derivatives depending on the skew angle. The purpose of 
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this analysis is only to localize those derivatives whose contribution to 
the root location is negligible and to compare them with each other. 

For the sake of clarity, we have labeled the roots according to the 
classical definitions though, in the case of a skewed wing, they may lose 
part of their meaning. Table 5-II reports the results of this investigation 
and shows the influence on each natural mode. They were evaluated according 
to the percent change in the natural modes while the value of the derivative 
ranged from zero to the maximum value assumed. The symbols used are: 

: less than 5 $ 

* : 5 to 20^ 

** : 20 to 50$ 

*** : 50 to 80$ 

**** : over 80^ 



In section 5.2 we have analyzed the spanwise lift distribution of the 
aerodynamic forces as computed with our numerical models. The stability 
derivatives are the direct consequences of those spanwise distributions . 


To date, neither program computes the stability derivatives with .. 
respect to rate of change of be . These derivatives take into account 
the time required for the effect of downwash produced by the wing to reach 
the horizontal tail. For symmetric aircraft, the methods described in 
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References 27 and 28 are usually used. 


The effect of downwash is to reduce the angle of attack of the 
stabilizer. Because the wing is at a skew angle, the distance from the 
wing to a conventional stabilizer varies from tip to tip and complicates 
the computation of such derivatives. 

Since the analysis done in section 5.3 showed that the influence of 
these derivatives on the dynamic of the aircraft is negligible and because 
of lack of time, we decided not to further investigate the analytic and 
numerical evaluation of such derivatives. 

We recall that the nondimens ionalizing quantities used in our 
analysis [Appendix A] are: cruise speed, wing span, and mean aerodynamic 
chord, both for the unskewed configuration. 

The results obtained, together with the Boeing ones, are shown in 
Table 5-III. Tables 5 -IV through 5-IX are the computer printouts of the 
first six cases of table 5 tIII. The symbols in these printouts are 
explained in Appendix E. 

In addition to the stability derivatives, the wing contributions 
are reported separately. 

Whenever a stability derivative differs from its corresponding wing 
contribution, this is because of the tail contribution. 

We shall now make a few remarks concerning the most influential 
derivatives (Table 5-II) as computed with lifting line theory. 

The comparison between rigid and flexible wing, in terms of stability 
derivatives, can be done by comparing columns 3 and 4 in Table 5 -III. 
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TABLE 5*11 - INFLUENCE OF THE STABILITY DERIVATIVES ON THE NATURAL MODES. 
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* No Stability Augmentation (SAS) included. 


** These data were taken from Boeing study [Ref. 5] since our program 
does not have the capability, to date, of computing these derivatives 


The eight columns represent: 

1) Strip Theory: spanwise flat and total lift distributions as from 
Reference 5 (Fig. 5.35). Sideslip evaluated according to empiri- 
cal method described in section 2.5. 

2) Strip Theory: spanwise flat and total lift distributions as com- 
puted by lifting line theory (Fig. 5.2). Sideslip evaluated as 
for case 1 . 

3) Lifting line theory, rigid wing. 

4) Lifting line theory, elastic wing (E.A. coinciding with quarter- 
chord axis. 

5) Same as No, 4, but no built-in twist. 

6) Boeing results rigid wing. 

7) Boeing results elastic wing. 
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The elastic wing, unlike the rigid one, shows a positive C 

mQi 

or pitch-up moment. 

There are no major discrepancies between our results and 

Boeing’s except in the derivative. Since the major contribution 

to the latter is from the fin, the difference may result from our 

a6 

assumption ^ = 0. 

It is useful to consider stability derivatives as functions 
of skew angle before discussing dynamic response. 

The effect of skew is seen in the accompanying graphs (Figures 
5.35 through 5.41) on a configuration that was stabilized by moving 
the C. G. .15 M.A. C. further forward. 

The behavior of these derivatives is explainable in terms of 
two important changes in the lift distribution resulting from an 
angle of attack perturbation. 

The first of these is the predominance of lift on the trailing 

wing. The second is the effective reduction of the lift forces 

arising from roll and pitch rates or changes in pitch attitude, that 

is, C T „, . 

1A is reduced. 

The program approximates as the calculated lift divided by 

dynamic pressure, wing area, and cruise angle of attack. From simple 

two dimensional sweep theory, the lift on a wing at constant angle 

of attack and free stream velocity varies as the cube of the cosine 

of the sweep angle. The effective dynamic pressure is reduced by a 
o 

factor Cos A and the effective section angle of attack by a further 
factor of CosA . The lift distrubition calculated by the program 
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closely matches this Cos A variation. C is proportional to C 

La L 

3 . 

and hence behaves as Cos A as well. 

It is not possible to easily follow the details of a change 

of skew angle through the lifting line analysis to verify the 
3 * 

Cos A dependence which results. 

The trailing vortices move closer to the downwash control 
points, and bound vortex segments decrease in length in proportion 
to CosA , which accounts for a factor Cos A . The change in the 
geometry relating the vortex segments to the control points is 
complex however and gives rise to three dimensional effects such 
as the lop-sided lift distribution, which are totally unaccounted 
for in two dimensional simple sweep theory. 

In addition to decreasing C , the closer lateral proximity 
of the skewed wing to the body reduces the roll induced velocities. 
However the extension of the skewed wing fore and aft of the C. G. 
causes increased pitch induced velocities. 

In pitch, roll, and yaw manouvers the lift perturbation 
distribution is unsymmetrical with respect to the center line so 
that consideration must be given to the changing moment arms in 
interpreting the behavior of the derivatives. 

The moment arms are decreased for rolling moment derivatives, 
increased for pitching moment derivatives and relatively unchanged 
for yawing moment derivatives. 

These observations lead to the following understandings of the 
effect of skew on individual derivatives: 

(1) C^, C T/V experience a decrease as explained earlier. A 

3 

Cos A dependence applies only to small variations about 
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an equilibrium position. The airplane would be maintained at 

a constant cruise velocity by increasing cruise angle of attack 
-3 A 

as Cos A. . 

(2) : The tendency of lift to shift to the rear wing leads 

to a rearwand shift of the neutral point as shown in the graph 

of vehicle aerodynamic center, or an increase in the static 

margin. However the decrease in C,„ dominates, and C = C x 

m mot m 

static margin is reduced. 

(3) : displays slight increase until skew reaches 20° decreasing 

thereafter From Appendix B, C _ = C (l-2kc ) C T and C_„ are 

xoc l m l La 

both decreasing which results in the calculated reversal. 
("Reversal" will be used to mean a reversal in the sign of a 
curve’s slope, and "reversal in sign" a zero crossing of the 
curve). The second term in the parenthesis is the contribution 
of induced drag, while the first results from the rotation in 


the perturbed body frame of a lift vector fixed in stability 
reference axes. 

(4) C : increases with the appearance of a sideways rotation of 
the lift vector and a sideways component of induced drag. The 
reversal is caused by the dramatic decrease of C and C at 

Li LCt 

high skew angles. 


(5) C = -C 

za. . ux 

(6) C^, where & is rolling moment, increases with the shift of 


lift to the trailing wing, but eventually reverses with 


decreasing and moment arm 

(7) C is the result of X , the X force due to an 0L perturbation, 
Tp. 0L 

acting through a moment arm proportional to CosA and through 
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an arm proportional to Sin A : initially side force generated 


yawing moment predominates and C increases. Decline sets 

TjCX 

in as the longitudinal force generated moment becomes in- 
creasingly smaller and C itself begins to decline. Again 

the reduction of C, and C „ with skew has an important role. 

L LQ; 

(8) C : Monotonicly increases for the same reasons as C . 

yp ya 

The increase is less dramatic because of the decreasing distance 
of the wing tip from the centerline, and reduced roll induced 
angles of attack. The non-zero value at aero skew is due to 
the rudder. 

(9) C : Normally zero for symmetric aircraft which have an 

zp 

antisymmetric roll generated lift distribution, the C ^ 

curve has two startling reversals and a reversal in sign. 

The derivative is very sensitive to slight departures of the 

lift distribution from antisymmetry. In the range 0° to 10° 

skew, lift decreases more on the leading wing in a roll than 

it increases on the trailing wing. After 10° the trend 

is reversed until the lift increment on the right trailing 

. o 

wing dominates as in figure 5-18 for A = 45 . 

The final reversal reflects decreasing lift curve slope 



curve slope cause a reversal. 
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(!2) C^ p : Decreasing X- force yawing moment anus are offset by 

increasing Y- force arms. The longitudinal position of the 

center of gravity will influence the shape. Figures 5-19 (a) 

and (b) show distortion of the expected antisymmetric X and Y 

force distribution. Since the incremental lift distribution 

in roll shown in figure 5-18 for A = 45° is almost normal, the 

distortion is related to an abnormal downwash distribution, 

figure 5.13 a). The departure from an antisymmetric force 

distribution decreases C from 0° to 30°. 

IP 

(13) C : Figures 5-22 (a) and (b), showing pitch generated 
X and Y- force distributions, are similar to those for roll 
generated horizontal force distributions. figure 5-21 the 
corresponding lift increment is noticeable greater on the 
trailing wing. Skewing the wing has a less detractive 
effect on pitching stability derivatives because induced 


velocities and aerodynamic twist increase. The detractive 

effects of diminishing effective angle of attack, and dynamic 

pressure remain however. The horizontal force distributions 
o 

are symmetric at 0 skew but become less so as portions of 

the leading wing extend in front of the center of gravity. 

The distortion mentioned earlier opposes this tendency. 

(14) C : declines steadily with an increasingly antisymmetric 
zq 

lift distribution (Figure 5-21) 


(15) 


C» : increases for the same reasons that 


C decreases, 
zq 


until the familiar high skew angle deterioration takes place. 
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(IS) C : is initially non-zero for a center of gravity not on the 

quarter chord. The introduction of antisymmetry in the incremental 
pitching lift distribution, and increased induced velocities, and 
pitching moment arms increase C rapidly at first. 

(17) C : T ^ e rap; *-^ increase after A = 15° , and the tapering 

off above A = 30° correspond to the varying antisymmetry 
in C . The initial decrease may be because the wing quarter 

q 

chord lies behind the center of gravity at A . 0 , so that 

when side force components of induced drag first appear their 

contribution to yawing moment is positive. As the skew 

angle increases, and the forward wing extends ahead of the 

center of gravity, its contribution to yawing moment 

becomes negative. The non-zero initial value of C could 

nq 

not be explained, but may result from an unsymmetrieal 
downwash distribution at zero sweep. 
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Figure 5.35 - Lift coefficient, lift curve 
slope, and aerodynamic center 
versus skew angle 
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Figure 5.41 - u stability derivatives 
versus skew angle 




5 . 5 Dynamic Results 


5.5.1 Natural Modes 

In 5,1 , Model description , it was noted that a C. G. location 
of 355 M.A.C. resulted in longitudinal instability. In case 3 of Table 
5-III the instability is not apparent from which is less than zero, 
however the roots illustrate the unstable short period mode discussed 
in 5.1: 



Damping 

Wn 

Real Root 

Imag. Root 

Phugoid 

. 1265 

.063 

-.00797 

+ .0652 

Dutch Roll 

. 1001 

1.058 

1059 

± 1.053 

"Rolling-Short-Period" 

.8007 

2. 120 

-. 1713 

+ 1.25 


Spiral Mode .0166 

Unstable "Short Period" . 1542 


Table 5-X records similar roots for case 1 and case 5. These 
configurations are therefore longitudinally unstable as well. The stan- 
dard short period and rolling convergence modes are recovered for both 
cases when 40$ SAS is introduced (equivalent to moving the C. G. forward 

.4 M.A.C.) The added static margin decreased C to -2.85 and -1.35 for 

ma 

cases 3 and 5 respectively* 

Since the numerical integration is sensitive to the number 
of wing stations, to obtain good results a large number of stations 
should be considered. 

For example, the stability derivatives for cases 3 and 5 
(40$ SAS) are shown in Tables 5-XI and 5-XII and the corresponding- 
natural modes in Table 5-X for 36 wing stations, compared to 40 for 


previous cases. 
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The roots of the dynamic equations were computed at 


.000 o 

A = 0 , 15 , 30 , and 45 for the configuration having a C. G. at 
-.15 M.A.C. , and the accompanying root loci plotted. (Figures 5.42 
through 5.46) 

These reveal the variation of dynamic behavior corresponding 
to the stability derivative versus skew angle curves plotted earlier. 

The effect of skew is to cause minor variation in the 
natural frequencies of all but the rolling convergence mode, which 
experiences a large reduction in roll damping. This is traceable to 
the diminishing of and further to declining moment arms, and 

induced velocities. 

The mode shapes corresponding to the new roots at A = 45° 
are shown in the form of Argand diagrams. (Figures 5.47 through 5.50) 
Although the characteristic roots have changed only moderately the 
dynamic modes no longer resemble the familiar ones of a symmetric 
airplane . 

The diagrams show the phase and magnitude relation of the 
six state variables describing the dynamic state. When the magnitude 
of a state variable is insignificant its relative phase angle is still 
shown by a line. State variables for non-oscillatory modes have only 
0° or 180° phase relations between them. The rolling convergence mode 
is one of these. The Argad diagram for A = 45° shows the diminished 
roll rate r resulting from the reduction in roll damping. 

In summary, the most dramatic change in dynamic character- 
istics appears to be the reduction in the roll mode damping. No 
other large changes in root locations due to sweep were encountered 
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for this configuration. The changes in mode shapes due to the long- 
itudinal/lateral coupling are generally small with the exception of 
the short period; which takes on a substantial amount of rolling along 
with the normal 0 and Oi . 
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Figure 5.46 - Root locus versus skew angle for 
lateral oscillation or Dutch roll 
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Rolling Convergence Mode 


Spiral Mode 
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a°= -1.063x16 _ 
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u/u =-1.260x10 
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0L =-4.132x10 
8 = 1.299x10 
/S = 4.509x10 
r = .039 

<p = 1.0 


Figure 5.47 - Effect of skew on rolling convergence 
and spiral mode shapes 
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Figure 5.49 - Effect of skew on short 
period mode shape. 
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A = o° 


C p = .465 @ 8.3 
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A = 45 C 


Figure 5.50 - Effect of skew on lateral 
oscillation on Dutch roll 
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5.5.2 Transient Response 


The transient responses have been computed for two cases: 

(1) Rigid wing, 40^ SAS added. 

(2) Flexible wing, no twist, 40 $ SAS added 

Figure 5.51 compares with 0! and j3 responses to an aileron im- 
pulse for the 40 $ SAS rigid wing case and shows the cross coupling 
between the lateral and longitudinal modes. 

Figure 5.52 compares the 0C and /3 responses to an elevator im- 
pulse for the elastic wing with 40 °jo SAS. Note that the /3 response 
is larger than the a response after the initial peak in the first 
second and will no doubt have some effect on a pilot’s evaluation of 
the handling qualities. This same result was found for the rigid wing 
with 40 io SAS. 


ORIGINAL PAGE IS 
OF POOR QUALITY - 
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RIGID KING 40% SflS RESPONSE TO AILERON 
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5.6 Summary 


The model used arid the aerodynamic result obtained with both strip 
and linear theory were discussed in detail and compared. The strip 
theory approach, though very simple, showed results not always satisfac- 
tory; in addition, its limitation to rigid wing case only turned out to 
be too restrictive, since the elastic wing behavior strongly differs 
from the rigid one. 

The effect of stability derivative changes on the natural modes 

was tabulated. The derivatives of a longitudinally stable configuration 

o o 

were graphed for skew angles from 0 to 45 , and discussed in light of 
the geometric and aerodynamic effects of skew. 

Substitution of these in the eighth order dynamic system gave the 
characteristic root loci as functions of skew angle. 

Calculation of the mode shapes revealed significant changes in 
dynamics, although the natural frequencies had not changed greatly, 
(apart from the rolling convergence mode). 

Transient responses to elevator and aileron deflections were 
presented. 
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VI CONCLUSIONS 


The analysis and computer program for the adaptation of the lifting 

line aerodynamic theory to the oblique wing have been described. 

The existence of a side force component, due both to induced drag 

and to the tilting of the lift force when the wing is skewed, has been 

shown. Stability derivatives were obtained for an aircraft used in a 

o o o o 

Boeing study with the oblique wing placed at 0 , 15 , 30 , and 45 . 

The derivatives were generated using the lifting line theory and a 
simple strip theory (for sweep = 45°). The two results are compared 
with each other and with the Boeing results and show reasonable agree- 
ment in most cases 

The stability derivatives computed using the lifting line theory 

were used in a linearized dynamic model of the aircraft to determine 

the effect of sweep on dynamic behavior. No instabilities or large 

changes occurred in the root locations for sweep angles varying from 
o o 

0 to 45 with the exception of roll convergence. The damping of the 
rolling mode was reduced by more than an order of magnitude due in 



A dramatic increase in the characteristic roll angle, in comparison 

to other state variables, was prominent in the rolling convergence, and 

. o 

the three oscillatory modes at A = 45 . The rolling motion in the 
Dutch roll is exaggerated with increasing skew, and surprisingly both 
the phugoid and short period modes picked up significant rolling motion. 
In the latter mode rolling dominated by a factor of three . 
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APPENDIX A 


A. I Relations Between Dimensional and Nondimensional Derivatives (Sta - 
bility Axes ) . 

The geometric quantities b^ , wing span, and c , mean aerodynamic 
chord, are referred to the wing in the unswept position. 


x = 2 ^ Cx 

“ u o u 


lb • sec 
ft 


x a ‘ I s Cx a 


[lb] 


x a “ qs 2^ Cx a 


[lb * sec] 


Xg = - mg cos0 Q [lb] 


x p = qSCx p 


[lb] 


x r qs w 0 cx j 


[lb • sec] 


*P * " S 2^ % 


[lb • sec] 


x q ■ qs sr % 

u 


[lb • sec] 


K t - qs J5J- CK r 


[lb * sec] 

The same relationship can be obtained for y and z 


k. 
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L a “ qSbc 4a ; 


M a = qScCra a 


L & " qSb 2U Q C £Cx ; 


M a = qSc 2U^ 


L p ' qSbC j£p 


; Mp = qScCnip 


L p - qSb 2^ % '• M p * 1 Sc 2^ °”p 


L p = 1 Sb 2^ % ’ M p * " Sb W Q % 


L = qSb C 


2U q Xq * 


M = 


q SC ^ ~ % 


K = qSb ^r- C 


2U Q "Xr ’ 


- b 

M = qSc rrr- Cm 
r ^ 2 U q r 


The relationships for the yawing moment are 
rolling moment. 


[lb • ft] 

[lb . ft • 
[lb . ft] 
[lb • ft • 
[lb . ft • 
[lb • ft • 
[lb • ft • 
the same as 


sec] 


sec] 


sec] 


sec] 


sec] 


for the 
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APPENDIX B 


Since only the wing contribution to stability derivatives is a 
new element in the stability analysis, we shall limit ourselves to deriving 
only wing derivatives and provide a list of the ones which can be com- 
puted in a classical way and can be found in the literature [Ref. 15, 21 
and 22], The normalizing quantities b^ , c are referred to the A- 0 
condition. 


B.l "u" Derivatives [Ref. 22], 




_ 1 _ 

qs 




+ M 

a 



qS m 
n a 



(B-l) 



C 


z 


u 


_1_ 

qS 




ac_ 


+ M 


a m 


(B-2) 


B,2 "p" Derivatives. 

Figure B-l shows the case of the aircraft experiencing a sideslip 
velocity v . 

The effect of introducing a p corresponds to a - A A for the wing. 
The fuselage also is affected, but we neglect its contribution since it is 
too complicated to evaluate it. 
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Neglecting the induced drag due to the fin, we consider only the 


component of the lift due to the fin in the direction. 


= - D cosp + s i n (P " 5) 


r -I1 F = 

x D S dp x 

p 


3C S 

cosp + C D Sinp + — C L (P " 6) sin (p - 6) 

% 


'in 


+ C cos (p - 6) ( 1 - 
L Fin \ 


p=o 


Since when p=0 also 6=0 

dC 




^ +C U X 


Sp / sp 


bC bC. 

D ^ D 


3C 


D. 


dp ~ UK dA 


ac. 


d, 2c ac T 0 . . 

L L 2 _d_ _ 1 _ 

aA tt^r aA l aArTiR 


(B-3) 


since 




cosA 


(B-4) 


and assuming 


we obtain 


<C L>A ' (c l\ 


cos 2 A 

cosA^ 


'&C. 


L 

dA /A=Aq 


tt Ai 


tanAg 


(B-5) 


(B-6) 


* C = 0 assuming the trimmed condition does not require any rudder. 
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I t,n A) 


(B-7) 


\ 3A TT& 




and after substituting into (B-3), the final expression becomes 


be 3c 2C? 
c x " ' ap ' aA " " rtiR "A) 

p 

or, in terms of the induced drag coefficient 

c *„ ~ ■ 2 V UnA o 


(B-8) 


C 

Jl 


In a similar way we can compute C . The side force coefficient 

A 

for a sideslip is 


C = - C cos8 - C sinp 
y c K D K 


(B-9) 


where C c derives from fuselage and tail contribution 


ac ac 

- C c sinp + + C D cosp + 


3C 


P=o 


(B- 10) 


~ + C 

ap d 


The main constribution fo usually comes from the body and the 

vertical tail, for the oblique wing case the side force generated by the 
wing also should be included in C c , but it turns out that such contri- 
bution is of second order and therefore negligible. 

The tail contribution is a conventional one and it is given by 
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Z£ . c (i . M 
^ Fin Wn' * 


(B- 11) 


similarly to the C derivatives, no attempt is here made to evaluate 

X P 

the fuselage contribution. 


3c 3C 

SA “ ' 2 c l tanJ V> 


(B-12) 


For simplicity, we shall neglect the influence of fuselage on C, 


Wing Contribution 


(— ) 

' Sp /w 


(B- 13) 


The magnitude of this term is strongly dependent on the FWLD, there- 
fore no approximation is possible which would give a meaningful result. 
Assuming a zero rolling moment for the cruise condition, it is possible 


to evaluate 


(sj 

' R ' \ 


by computing the wing rolling moment at the per- 


turbed sweep angle and dividing it by the increment in sweep. 


Tails (Fin) Contribution [Ref. 21] 


l SB /F 


S 1 

c L (p - 6) cos 0 - 6) — ~~- 

a p 


(B-14) 
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-2 
sp /f 


b5 M ^ZF 
Sp S / b 


(B-15) 


The same approach used in computing jC 


evaluating the wing contribution to I C 


is necessary when 


and C 


The tail contribution is negligible, therefore 


c = c 

im \ m„ / 

p \ p/w 


(B-16) 


The fin contribution is given by [Ref, 21] 


where 


n p 'Fin L a'Fin 


S„. 

Fxn Fin 


(B-17) 


(B-18) 


vertical tail volume. 


The sidewash factor , generally speaking is difficult to estimate 
with engineering precision. Suitable wind-tunnel tests are required for 
this purpose. The contribution from the fuselage arises through its 
behavior as a lifting body when yawed. Associated with the side force 
that develops is a vortex wake which induces a lateral-flow-field at the 
tail. The contribution from the wing is associated with the asymmetric 


structure of the flow that develops when the airplane is yawed. This 
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phenomenon is especially pronounced with low-aspect-ra.tio swept wings. 
When such tests are not available. References 27 and 28 can be used for 
empirical values. 

B.3 "OS' Derivatives, 



where 


q s ~i- — ~ = c 

x a qs ba. L 


SC D 

ba C L ( 1 ‘ 2k °L 


a 


Sc 2 


(B- 19) 


a 


For a rigid wing, ^ = tanA , therefore integrating (3.28) over 
the span and normalizing the result 

b/2 


. 3F 

c = J X = 

y a qs ba qs 


k f ( L a y 0 + L o tanA ■ 2k L o c l ,) tanA 

J u /o ' 


we obtain 


C 

y a 


-b/2 


C L ^0 +C L 1 ' 2kC L 

a x a 


tanA 


(B-20) 
(B-2 1) 


C 

JS; 

In stability axes 



(B-22) 
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m 


a 


The presence of an unsymmetric lift distribution due to the angle of 
attack introduces a new wing contribution to C 


m a 

Let us first evaluate this contribution before considering all the 

other ones which are common to symmetric aircraft. We shall refer to 

Figures 3.5 and B-2 for the symbols used. To be consistent with the 

conventional notations we shall consider this contribution as a part of 

C = C (A) referred to the M.A.C. at zero sweep (y , ). 
m Q nig r W MAC 

The wing moment about y^^ is given by 


(M) 


w, 


b/2 

/ 

-b/2 


Z y MAC ^ 


(B-23) 


where (M) is the wing pitching moment due to the effects of skew. 


w 


A 

Let us now consider the general problem. 


Moments about C.G. 

a) Wing contribution 

M = Mq + (M)^ + (L cosQp + D sina) (h - h ) c (B-24) 

A ■ 

b) Tail contribution 

M t = - \ L t (B-25) 

c) Total moment 

= M + (M) + (L cosa + D sina) (h - h )c -A. L 

iuw* U wbtt 
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Figure B-2 Moment about the C.G. in the Plane of Symmetry. 

N.P. = Neutral Point 

c = M.A.C. at A= 0 

OIq = Angle of attack (wing-body) 

€ = Effective wing downwash at tail 

i = H-tail trim angle 

t 

a fc = a - (e - i fc ) = H-tail angle of attack. 
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M_ + (M) + L 

0 w . 

A 


cos a 4- “ sina 

U 


(h - h ,)c - i L 
v wb t t 
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- M„ + (M) + L (h - h ) c - it L 

0 v 'vij^ v wb t t 


cosa 0 + L Sina 0 ~ 1 + L a 0 1 


qSc 


- ^ 


Vt 


C = C 4- (C ) + C T (h - h . )- C T 

m T m Q m W A L wb L t g - 


= C + (c ) + C T (h - h , )- V__ C 7 

' m'w A L wb HI 


where 


v 

H Sc 




(a 0 + *t 


e) 





V. 


H 




(C ) 


m w. 


+ C 


t ( h - h J 

L a wb 


1 - 


a 


5s 

3a 


H 



(B-26) 

(B-27) 

(B-28) 

(B-29) 

(B-30) 

(B-31) 

(B-32) 
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Neutral Point (C.G, location about which C = 0) 

m 

a 


C = ~~ (C ) + c (h - h )- V c 

m a Sa mw A L a wb H L a 


l - 


5e 

3a 


= o 


h - h = h - ~ 
wb C 

J u 


a 


3a (C m' ) w^ “ v h C L a l 1 “ 3a 


or 


h „b- 1 + cT Ta (c Jv. - v h c l„ {' - Is 

L a 1 A a t 


and by substituting (B-34) into (B-32) 


C = ~ (C ) + C T (h - h ) 

"a 5a " "A L a 


act ” v h C L ( 1 " 3a 

A a v 


- v h c l q P'fs 

t 


C = c (h -h) 

m a L a 


h - h - Static Margin 


Aerodynamic Center Build-up. 


where 


X. =h c 
AC wb 


h = .25 

wb 


a 


X AC^ V H C T 

Tails L 


(l- h) 

\ 3a / 


a 


(B-33) 

(B-34) 

(B-35) 

(B-36) 

(B-37) 

(B-38) 

(B-39) 
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i a 


(B-40) 


X 


AC. 


Pivot and Sweep 


3a ^m^w 


a 


A 



This is one of the new derivatives and only the wing contributes to 
it. It is formed of two terms: one due to the lift and one due to the 
side force. 





(X x + Y y) dy 


(B-41) 


The knowledge of the spanwise distribution of X and Y is needed 
in order to evaluate C 

n 


B.4 "p" Derivatives. 

The local angle of attack varies linearly according to 




(B-42) 


For the symmetric wing case this derivative is zero since the in- 
crease in lift on one side of the wing is balanced by an equal decrease 
on the other side. 


This is no longer true for the oblique wing, therefore 



2U_ 


Sb . 



(B-43) 
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The evaluation of 


dZ 

dp 


can be done according to the method outlined 


in section 3.2.2. However, such an approach will require the evaluation 

ax 

of — that will take into account the behavior of the oblique wing. 

To remind us of this detail, we shall add the subscript p to the 
derivatives that must be evaluated in that way. 


Since 


and 


. s 5a 

JL _ Jl p 

dp da dp 


da - 

E. _ JL 

ap u Q 


dZ 

dp 


_ /dz\ x /&& \ 

\da/u " [\daj 


+ d. 


X 


(B-44) 


Therefore 


- 2 

Sb„ 


/ 

J -b [2 


dSj + d o 

p 


x dy 


(B-44) 


Wing Contribution. 
Similarly to (c 


was antisymmetric. 


, , ,C would be zero xf the change xn Ixft 

z \y 

p-'w ' ■'p'w 


2U 


0 


b/2 

/ * 2 U n 

^ dy = V 

-d/ 2 U '-b/2 


b/2 

L 


dp 


(1 y - d tanA) dy 


(B-46) 


180 



We notice that 


■^ = 0 for a rigid wing, but it would be different from 
zero for a flexible one. 


Fin Contribution [Ref. 21]. 


( Cy J: 


y p/Fln 


f!»c 

b ° s V 


U, 


£6 Jb_ 

3p 2U. 


& • s 

FZ F 


ll n l 7 ^. _ ^ 

s L a . \ b 0 Sp 

Fm 


Therefore, for a rigid wing 


(B-47) 


y p b o s 


b/2 

f [(fa) y o - (Is) tanA 
4/2 1 p p 


dy ' T c i 


<U 


cl . v 

Fin 


fz _ S 6 

5p , 




(B-48) 


For the flexible wing, the term ^ is no longer zero, therefore 
the term should be included in the integration representing the 


wing contribution. 


s 

— E 


Wing Contribution 



b o s 



sz - . 
tr xd y 


- 2 

b o s 




dy 

(B-49) 



Fin Contribution 



C 

m 

E 


No tail contribution to this new derivative 



(B-52 ) 


C 

n 


Wing Contribution 

The side force is the new element in the wing contribution to C 

n 

P 




Tail Contribution [Ref. 21], 


/ \ 

FZ Sp 

C = V 

\ n p ' Fin V 

L b 0 Sp J 


(B-55) 


°kn 


Therefore, by adding the two contributions, for a, rigid wing we 
obtain 


n 


b/2 

/. 


V "-b/2 


JL - (— 

*0 \$a 


x + 


+ v 


V 


21 


FZ S6 
>0 ” Sp 


I) hi 


iV 

Fin 


5d 

Sa 


tanA 


y j dy 

(B-56) 


The flexible wing would have the extra term 


l 

0 3a y 


B.5 "q" Derivatives. 

The q derivatives are derived in the same way as for the p ones. 

All the comments made in the previous. group of . derivatives can be 
extended to this one and therefore they will not be repeated. 

The local angle of attack now varies linearly according to 



Wing contribution only 


a 
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/ \ 2U„ 

C y “ (c ) ■ — 

y q V ^q'w cS 


b/2 

2 /* 

CS 4/2 


fca/q 


b/2 


2ij 

b/2 



/ 

$ dy 

= 

cS 

/ (»)i 


“'-b/2 


-b/2 







(B-57) 


b/2 


y 0 ' (fOq 



_2_ 

cS 

/ 

-b/2 

/di\ 

\^/P 

tan A 

y d y 


For flexible wing the term JL y must be included in the integra- 

U CM" "" 


tion. 


C 

z 

__a 


Wing Contribution. 


2_ 

cS 


b/2 

/ 

-b/2 


q w cS 


— ) y dy = 

aa/q 


b/2 

/ 

-b/2 

b/2 

CS 4/2 


dy 

aq y 


.M 

,a»/q 


(B-58) 


4- dO 


y d y 


Horizontal Tail Contribution [Ref. 21] 


Vt 


and 


J q/H.Tail Sc 

b/2 

In 


V = " 2v h c l 


a 


HT 


cS 


M , , H 

» q + d 0 


y dy - 2 V H C L 


(B-59) 


(B-60) 


^HT 
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Wing contribution only 


V 2U_ 

b/2 

/• 

C - ° / 

^q'w b^ cS 

-b/2 

b/2 r , 

/ U 

+ d„ 

J [W<1 

0 

-b/2 



to 

to 


x dy 


x y dy 


m 


Wing Contribution 


( c „) ' 
V a /w 


q ' w c S 
b/2 


2U 0 »- H 
^ si ydy 


-2 

c 2 S 




1 2p;/q 0 


'-b/2 

Tail Contribution [Ref. 21] 


-2 ^ 

y dy 


\ z t 

C m =- V H C L “ 

q/T c 


Therefore 


m 


b/2 

-2 r 

- 2 . / 

C S -b/2 


4. A 

toe) q o 


-2 ^ 
y dy 


V H V 


a c 
HT 


(B-61) 


(b-62) 


(B-63) 
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c 

n 

__a 


Wing Contribution only 



(B-64) 
y dy 


B, 6 "r" Derivatives . 

Except for the new C derivative, due to the wing aerodynamic 

r 

coupling, and for the wing side force contribution to yaw moments, all 
the other derivatices are standard. 

The method used in evaluating r derivatives in section 3.2.4 will 
be applied here without any further explanation. 

Assuming the steady state flight condition to be straight levelled 
flight 

c = C = C = 0 
z y x 

T* r* T.f -r* 


C 

JjE 

Fin Contribution only [Ref. 21] 




2 L 


FH 




m 


(B-65) 
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Wing Contribution 


_9TT h > 2 b/2 

(' h ) ‘ rrf- / h * ^ dy = ^ / 

\ r/w qb Q S J h/Z b n S 


b 0 S *'-b/2 


—2 

■ 2 q x 


Tail Contribution [Ref. 21] 


'Z. 


= C 


r/Fin 


^FZ 
b 0 


Therefore 


-d 

V 


b/2 


-2 


FZ 


X n x dy + C , 

0 V b n 

■'r 0 


m 


Wing Contribution only 


C = C 

m r V m r /w q b 


9n b/2 b/2 

) / 

/w qb n cS y v/0 01 b Q cS J hj2 


0 -b/2 


x Q yx 


Wing Contribution (including side force) 


n 


2U 


0 


r^w 


q b 0 S -b/2 


b/2 

f. 


d Q x + (^ 0 y Q - d tan/) y 


M I 

Sr 


dy 


(B-66) 
(31 - 67) 

(B-68) 


dy 

(B-69) 



b o s 


b/2 

/ 

-b/2 


d Q (x - tanA y) + JL 


0^0 


y) 


x dy 


(B-70) 


Tail Contribution [Ref. 21] 



(B-72) 
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APPENDIX C 


C. 1 Similarity Transformation between Inertia Matrices . 


The conversion of an Inertia Matrix from body into stability axes 
occurs according to a similarity transformation [Ref. 26]. 


where 


tv = £T 8/b ] (V n: a/b y 


<C-1) 


[I ] = Inertia Matrix in stability axes 

[I ] = Inertia Matrix in body axes 




[ ] 


= Rotation (or direction cosine) matrix from body to 
stability axes 
= Transposed matrix 


The stability axes are obtained by rotating the body axes about the 
y^ axes by an angle a , therefore y^ = y g . 

The rotation matrix [T g ^] i- s given by 


and 


[T s/b ] 

[X ,, ] T - 

L S/D J 


cosa 0 
0 1 

sinQ! 0 

” cosa 0 
0 1 

_ s ina 0 


sina 

0 

cosa. 

- sina 
0 

cosa. 


(C-2) 


(C-3) 


In a symmetrical aircraft both the I and I components of 

xy yz 

the inertia matrix are zero. When the wing is skewed, the I compo- 
nent is no longer zero, therefore for an oblique wing aircraft the 
inertia matrix, in body axes, looks like 



CV “ 


XX 


xy 


xz 


xy 

[ 

yy 

0 


xz 

0 

[ 

yy J 


(C-4) 


I£ we know transform [1^] according to the similarity transfor- 
mation (C-l) , we shall obtain the inertia matrix expressed in terms of 
the stability axes, [I g ] . 

After some algebra we obtain 


1 

I' 

I' 

XX 

xy 

XZ 

? 

I' 

I' 

xy 

yy 

yz 

» 

i' 

1’ 

'xz 

yz 

zz 


where 





2 


. 2 

I 1 


I 

cos a + 

I sin2a + 

I sm a 

XX 


XX 


xz 

zz 

I' 

— 

I 




yy 


yy 




i' 

= 

i 

sin a - 

I sin2a + 

2 

I cos a 

zz 


XX 


xz 

zz 

I’ 

— 

I 

cosQI 



xy 


xy 




I' 


I 

cos2a + 

x a - i 

) sin2Q! 

xz 


xz 


l ZZ XX 

I' 


- I 

sina 



yz 


xy 
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APPENDIX D 


D . 1 Downwa sh Matrix fS ^ l ■ 

The lift or circulation distribution can be visualized as resulting 
from a system of horseshoe vortices, each of which is of constant 
strength (Fig. D-l). 




Figure D-l [Ref. 14]. 
Horseshoe Vortex System 
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The net strength of the trailing vortex at any point on the span of 
the wing is numerically equal to the rate of change of strength of the 
bound vortex in the spanwise direction. Little loss in accuracy with 
respect to the spanwise airload distribution will be entailed if: 

a) The total strength of the chordwise system of bound vortices 
is concentrated in one bound vortex located at the local span- 
wise quarter-chord point. 

b) The downwash angle at each vortex station across the span of the 
wing, at the local streamwise three-quarter- chord point (= down- 
wash control point D ) is equal to the geometric angle of attack 
for airfoil having a 2-D lift curve slope equal to 2TT , When 
the section 2-D lift slope is different from 2lT , equation D-l 
must be used 

/W\ m o 

(v) 3C/4 _ 27T a f (D_1 ' > 

The downwash angle at any one downwash control point is the sum of 
the incremental downwash angles due to the horseshoes in the system of 
horseshoes which represent the wing and its lift distribution. 

Assuming the geometry of the wing platform, the angle of 
attack and section 2-D lift curve slope variation are given across 
the span, the unknowns arc the values of the running lift at each 
point on the span. The strength of each bound vortex represents 
the average airload over its own portion of the wing span. 

The method for determining the downwash matrix for an oblique 
wing is now illustrated. 
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In Fig. D-2 a system of horseshoe vortices and the associated down- 


wash control points are shown. 



°°* 



Figure D-2 


At the yj station the section lift curve slope will be iru and 
the angle of attack of the section zero lift line is a : . 

Since a linear relationship exists between the strength I 1 of a 
particular horseshoe vortex j and the downwash velocity W„ at the 
point i 


where K is a constant, 
vortex system is 


w. . = k. . r. 

ij ij j 


and the downwash at 


(D-2) 

i due to the entire 


w. = X). w. . = y.. k. . r. 

i *-• 3 xj xj j 


(D-3) 
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From equation (D-l), which expresses the relationship which must 

W 

exist between the downwash angle — at each control point, the wing 
angle of attack and the section lift curve slope for the wing 

station at the control point, the following series of equations result 


V 21T a i 


W, 


_2 

W 


“l 

2 IT a 2 


and in general 


W. 


m i 

2 IT \ 


and substituting (D-3) into (D-4) 


(D-4) 



V lj J 


(D-5) 


.th 

l 


The relation between the running load 
station is 


1 and the circulation at the 


b ■ 0*1 


(D-6) 


Equation (D-5) can therefore be rewritten in terms of the running 
load 1 


W 


i 1 


„ „ Y.. K. . r. = ~~r T., K. . 1. 

V V 1J 3 pV 2 "j 13 J 


(D-7) 


Substituting (D-7) into (D-4) results in 


m. 


PV 


1 • 

K.. 1. = — a. 
l] J 27T i 


or 


qm. 
r 


V. K. . 1. = — a. 

J TT i 


(D-8) 


(D-9) 
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And in matrix form 


[K] U = 


“^0 
ff ' J 


a_ 


where 


ir v . 


is a diagonal matrix. 


Defining a new matrix 


[S x ] - [K] 


Equation (D-10) can be rewritten as 


(D- 10) 


(D- 11) 


[S 1 ] \&\ = ja f | (D-12) 

The elements of the [S^] matrix are to be influence coefficients 
relating the incremental downwash angle at each control point to the 
intensity of the running lift over each increment of the semispan of the 
wing. 
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D.2 Evaluation of K. . Elements. 

• ij- 

The velocity induced by a vortex of strength F at a point 
be written as [Ref. 24] 


W 


P 


F (cosa - cosG) 
47TR 


P can 


(D- 13) 


where Oi and p are the angles between the direction of the vortex 
segment and lines joining the ends of the segment to the point as shown 
in Figure D-3. 



Figure D.3 - Finite Segment of a Straight Vortex Filament [Ref. 14], 

A plan view of the geometry of a typical horseshoe vortex is given 
in Figure D-4. 

In order to evaluate the incremental downwash velocity induced by a 
single horseshoe vortex it is convenient to consider the following three 
cases : 

1) Control point to the left of the horseshoe vortex. 
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2) Control point within the horseshoe vortex. 


3) Control point to the right of horseshoe vortex. 

Fig. D-4 shows the quantity that will be used in the derivation; by 
the subscript i we shall indicate quantities relative to the control 
point D. , whereas by j we shall refer to the horseshoe reference 
point Vj . 

The origin, of the axis system is at root quarter chord point. When- 
ever the locus of the quarter chord point does not lie on a straight 
line, it is assumed to be given by an equation £(y) w.r.t. a straight 
line passing through the root quarter chord point and aligned with the 
wing span (unswept case) 



X. 

2 


X. 

x 


C. 

i 


1) 


, Y. coordinates of horseshoe reference point V. 

2 *2 

C. 

- — , Y^ coordinates of control point 
chord length at station i 

X. = - y . tan A+ f(y t ) tan A 

Control point to the right of the horseshoe vortex (Fig. D-4a) 
defining 
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R. . = Ax. . + Ci/2 
ij ij 


AY ij - V Y j 


Fig. D-4a Control Point to the Right of Horseshoe Vortex. 
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For the left trailing vortex 


a = 0 


p = a + 90 


and equation (D-13) becomes 


[1 - cos(a. . + 90 )] 

(W ij^L = r j Ai(AY^Vh) 


1 + s in (X, 


JJ_ r 

» N X - 


4lT(AY i j+h) j 


(D- 14 ) 


sin d. . = 


R. . 


V (Ay. .+ h) 2 + R? 4 


IJ 


1J 


(D- 15 ) 


For the right trailing vortex 


a = 90 - p . . 

1J 


p = 180 


[cos(90 - p ) - cos 180 ] 
(W ij } R = _r j 4 tt(A Y„ - h) 


sin p + 1 

„ :! l .. lj 

4ir(A Y. -h) 


T. (D-16) 


The minus sign derives from having assumed positive downwash veloci- 
ties. 


sin p . . = 


R. . 


lj V (A Y. . - h) 2 + R 2 
1J 


(D-1.7) 


For the bound vortex 


a - a. 


ij 


p = B. . 
ij 

(cos CL. . - COS P . . ) 

w. . = r. 7^ U ~ 

ij j 4 rr r ± . 


(D-18) 
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Since 


i s j 

Ax.. = Ay.. = o 

X J 11 


R. . = C./2 
ii 1 


Both trailing vortices will induce the same downwash velocity 


a = o 


s = a. . + 90 
1 n 


< W ii>L = (W ii> R = < 1+si " “-<> 


11 


(D-23) 


sin a. . = 
n 


R. . 

li 


Ch/2 


V h 2 + R^' V h2 + (Ci/2)‘ 


(D-24) 


For the bound vortex 


a - a 


ii 


B — 180 + 0! . 
r ii 

r. 

(W. . ) = T'~a f COS a - • + cos a . . 

11 B 4ff \ 11 ii 


(D-25) 


Therefore 


cos a. . 

ii 


ypr-T-? 


h + ~r C . 
4 i 


(D-26) 


W. . 

11 


4tt 


2 4 

— (1+sin a..) + — cos a . . 
h x ii' C. ii 

i 


(D-2 7) 


3) Control point to the left of horseshoe vortex (Fig. D-4c), 
defining 
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For the left trailing vortex 

a = o 

p = 90°+ a: . 

ij 

rj.o n 



(D-28) 


<Vt = 


1 - cos (90 + a[.) 

AY. . , 

ij " h 


r, (1+ sin Ol! .) 
.1 13 

Ait (Ay„ - b) 


sin Ci\ . = 


ii 






ij 


For the bound vortex 


a = a. . = 

ij 


180 - a! . 

ij 




180 - B! . 

H iJ 


1 

(W. . )_ = t — (cos CC. . - cos B , . ) 

v ij'B ij p ij' 


cos oc ^ = cos (180 - ) 


Therefore 


AY. . - h 


cos ct. , = - cos a! . = - 


LJ 30 V (AY. . - h)^ + R?/ 

1J IJ 


S imi lar ly 


Ay. . + h 
JLL 


cos ‘ ' cos = ' “jw.w + a?.' 


ij 


ij 


For the right trailing vortex 


a = 90 - pi 


8 = 180 


(D-29) 


(D- 30) 


(D-31) 


(D- 32) 
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( Vr = 4tr 


r. [cos (90° - p ' . . ) - cos 180° ] r. (sin p ! . + 1) 

— . — 1 ij — ~ 1 . — — — 13 ■ — /n_-3 

4t T (AY^+h) 4 it (AY_+h) ii} 


sin = 


(AY. . + h) + R. . 

ij U 


Therefore 


r. 

w = — i 

ij 4 tt 


1+ sin a! . ) 1+ sin 8 . cos a. . - cos 8 . , 

- 12- . ill . 11 

(AY. . - h) Ay. . + h R. . 

13 ij ij 


(D-34) 


It is now possible to derive the matrix K of equation (D-10) 
since each element will be given by 


K = 

ij r 



D.3 Structures Fundamentals. 


As done in the previous section, the continuously varying spanwise 
airload distribution will be replaced by a series of constant intensity 
running loads . 

The assumption of the section aerodynamic center acting at the 
quarter chord becomes weak when considering the bending and 
torsion due to the airload distribution whereas it is quite good when 
computing spanwise lift distribution. This inconvenience ,due to the 
lack of predicting chord-wise lift distribution, can be reduced by 
introducing a correction factor f which will allow the section aero- 
dynamic center to be placed in any desired place along the chord. Such 
a correction factor is here assumed to be known. 


Because of reasons which will appear more evident in this section, 
the horseshoe vortices must be chosen in a way such that the aircraft 
centerline will coincide with one trailing vortex. 

Let us consider the geometry of the structural skeleton of the wing 
as shown in Figure D-5 and define 

M , M , ... , M Rolling moment due to total lift of all the 

X 1 x 2 X i 

vortices outboard of this point (positive when 
raises left wing tip). 

M , M , ... , M Pitching moment due to total lift of all the 

y l y 2 y i 

vortices outboard of this point (positive when 
nose up). 
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, ... , Beam bending moment at elastic axis point about 

axis perpendicular to local elastic axis (posi- 
tive when it compresses wing upper surface). 

, ... , T^ Torsional moment around elastic axis (positive 

when leading edge up). 


, L 2 , ... , Total lift acting on the wing section having 

span 2 h , numbered from the left wing tip 

L. = 2 h Si. 

1 1 

i being the intensity of the running lift at 
station i . 


6 1 ’ 6 2 ’ » e i Streamwise distance from horseshoe reference 

point at a wing station to the corresponding 
point on the elastic axis (positive when elastic 
axis point is to rear of horseshoe reference 
point) . 


f 1 , > * • • > ff Section aerodynamic center correction factor; 

(f positive when aircraft is to rear of quarter 
chord point). 


The general form for the bending moment is 

M. - M cos A + M sin A 
ix. y. 

1 J x 

and for the torsional moment 

T. = M cos A - M sin A 
1 y. x. 


(D-35 ) 


(D-36) 


At station 1, on the center line of the horseshoe vortex nearest to 



the left wing tip (Fig. D-6) the following equations apply: 

m * h. IjL = t 
M x x 2 2 L 1 4 


w _ 1 / , h tan A \ 

M y“T( e l‘ £ l °1 + 2 J 

and substituting into equations (D-35) and (D-36) 


M 1 2 


h cos A . 

+ (e 1 - f^ + - tan A) sin A 


(D-37) 


T i = T 


At station 2 


(e^ - f^ j tan A) cos A 


L h 

L i (h i +h) + — 


sin A 


(D-38) 


(D-39) 


= L 1 (AX 12 - t 1 C 1 + e 2 ) + L 2 (e 2 - ^ C 2 + - tan A) 


and in general, for the left wing (forward), for i>l 

i-1 

(“*.), -£ 2 h 1^(1 - k) + L | 

' i'L k=l 

i- 1 

(\)iT S + e i - *k V + f ( e i - £ i c i + 1 1 “ A ) 


(D-40) 


(D-41) 


(D-42)' 


where 


“kr^- x i 


herefore, for i>l , the bending moment and torsion for the left wing 
are given by: 
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Leading Edge 




Figure D-6 - Plan View of Left Wing Tip Section. 



l i ) L = Zj 2hL k (l-k) + ^ cos A 

1 I k= 1 


+ 2 L k ( ^ki + e i ' f k C k^ + ~2 (e i " f i C i + 2 tan ^ sinA 

_k=l 

'i-1 L 

Vl" 5 + e i- f k v + ~£ (e i- f i c i + f tanA) cosA 

re=l 


2 h " k) + ^ sin ^ 


For i = 1 expressions (D-37) and (D-38) should be used. 

For the right wing, assuming n to be the farthest right station 

L h 

M x = ff < D “ 45) 

n 

L h 

M = -£• (e - f C — 1 tan A) (D-46) 

y 2 n n n 2 
n 

h L h 

(M n ) R = - L n T COS A + "f (e n " f n C n " ~t tanA ) sinA 0>-*7) 


(T n } R = f (£ n ' f n °n “ ~2 tan A) COS A + f sin A 


and in general, for the right wing, for k<n 


It 1 

( M 4 = S 2h ^ (i - k> - L ^ 


(D-48) 


(D-49) 


L 

(vl - 'f + e i ■ f k v + t <e i - 4 h - ! c “ A) 

' i ' K k.-n 


(D-5 0) 
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where 


- \ ■ x i 


Therefore the bending moment and torsion, for the right wing and 


for k < n , are given by 


- E 


(M i ) = As 2hyi-k) - L ± ^ cos A 

Lk=n J 

i+1 ^ 

+ £ V®ki + e i • f k V + f (e i • £ i c i • I tan A) sin 

Lk=n 


( 0 - 51 ) 


i ) R * ? + \ ‘ f k V + T <e i ' £ i c i - 1 t,m A) 003 A 


Ti 2 h L^(i - k) - jr sin A 


(D-52 ) 


At the root station of the elastic axis the bending moments of the 
right and left wing cancel each other in a trimmed flight condition. 

The presence of an unsymmetric airload would introduce a non-zero aero- 
dynamic moment at the root station. 


In Figure D-7, such a case is shown. 

<Vl 



<Vr 


Figure D-7 
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where 


AM « (M r ) R + (M r ) L 

(M ) = Root bending moment due to right wing airload 

x R 

(M ) T = Root bending moment due to left wing airload 
X Li 

r = Root station 

In such a case the following considerations can be made; 

a) Instantaneously the wing is assumed clamped at the centerline, 
allowing for the discontinuity in the bending moment. The 
structural rotation angle across the span, due to bending, can 
be computed as shown in the next pages by simply substituting, 
for the semi-wing considered, the corresponding root bending 
moment . 

b) The AM produces an angular acceleration which can be decomposed 
w.r.t. the roll and pitch axes. 

This acceleration introduces angular velocities about the roll and 
pitch axes which results in an apparent or aerodynamic twist. This new 
twist contribution modifies the spanwise airload distribution and, con- 
sequently, the root AM. A change in AM changes the angular accelera- 
tion and so on. 

Because of this discontinuity at the root it is convenient not to 
consider the root station. 

In case such a point were included in order to take into account 
the two different values of the bending moment it would be necessary to 
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write two separate matrices: one for the left and one for the right wing. 
In addition, when computing the streamwise twist due to wing flexibility, 
this discontinuity in the root moment comes up again. 

By using a system of vortices such that the aircraft centerline 
coincides with one trailing vortex, the root double-valued point will 
not be computed. In doing so nothing is lost in accuracy but there is 
a gain in simplicity for the matrix analysis. 

In matrix form, equations D-37, D-38, D-43, D-44, D-45, D-46, D-51, 
and D-52 can be opportunely combined and become 

j M | = cos A [r] + sinA[u] jij (D-53) 

j T | = - sinA[r] + cos A[u] |l [ (D-54) 

where 



, JL are the values at the last station to the left of the air- 
craft centerline (Fig. D-8). 

, & are the values at the first station to the right of the 
aircraft centerline (Fig. D-8). 
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■P-IP* 







[A] 

! [0] 

[0] 

i l_n 
1 i -J 


where [A] and [B] are reported in Tables D-I and D-II. 

The moment at the aircraft centerline can be computed as follows 


A. ' \ + h 5 


M I = M + L (e - f C + - tan A) 


yJh L r LL 2 




K) 


- M + f (e r - f R C R - | tan A) 
K K. 


M ] = Rolling moment at the aircraft centerline elastic 

' r/L 

axis due to left wing airload distribution. 


M j = Pitching moment at the aircraft centerline elastic 
^r/L 

axis due to left wing airload distribution, 
rly (m^ ^ and j indicate the corresponding moment due 


Similarly |M ] and 
\ X r/R 


to the right wing contribution. 

The total rolling and pitching moments acting on the airplane are: 
for the rolling moment S. 


\ + \ + 4 <VV 
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£<v 


fjC^ + — tan A) 


( ^12 + e 2-W 


|(e 2 - f 2 C 2 +| tan A) 


0 




f l<V 


(Zsc 2 i + e i-.£ 2 c 2 ) 


7 (e l " f l C i + 2 tan 


O O 

*3 {33 

% g 

80 £ 

Sn 

J> >• 
t- Q 

£3 ^ 


^!L +e L- f l C l> 


( ^2 L +e L- f 2 C 2 > 


2 (e L' f L C L + 2 tanA) 


TABLE D-I 




TABLE D-II 



for the pitching moment M 


+M + i er ( L +V + | (^-ytanA- L tAVV R C 
7 R 

(D-62) 

The streamwise angle of attack contribution due to wing flexibility 
can then be obtained from 


where 


a 

s . 

X 



sin A + 



cos A 


(D-63) 


0! = Streamwise angle of attack contribution at station i 

S i 

due to values of bending and torsional moments acting 
inboard and at station i . 
k = R Station when i on the right wing 

k = L Station when i on the left wing 

El Effective beam bending stiffness around elastic axis 

GJ Effective torsional stiffness around elastic axis. 


Since dS = 


2h 


cos A 

for station 1 can be written as 


for this study, the integrals in equation (D-62) 


/ 


_M_ 

El 


ds = 




M 

h , 2 2h , 

. _ ~ 


E 1 I 1 cos A E 2 X 2 cos ^ 


M (L - 1) 2h \ 2h 
E (L-1) I (L-1) cosA e l x l cosA 


2h 


M, 


+ 


cosA 2 E-I., . „ 

1 1 x=2 l i 


~ e • i . 



I 

/ 


X ds ~ 2h 
GJ cosA 


2 Vi 


L 

+ £ 


1=2 Vi 


And similarly it is possible to obtain the expressions for the 
corresponding values at each station. 

In matrix form: 


where 


i a i.iL 
I si cosA 



[El] | M ! + cos A [GJ] j T } 


[El] = 


r cc] 

; to]' 

i — i 

S ! 

i — i 1 

i 

! ED], 


(D-64) 


(D-65) 


where [C] and [d] are reported in table D-III. 

The matrix [GJ] is similar to [El] and can be ontained from equation 
(D-65) by simply replacing El with GJ . 

Substituting equations (D-53) and (D-54) into equation (D-64) 


■ a l = XX 
s 1 cosA 


sinA [El] (cosA [r] + sinA [u]) 


(D-66) 


+ cosA[GJ] (sinA[r] - cosA[u]) 


which can be written as 


I a s I = 4h 


sinA([EI]+ [GJ])[r] + (^~~£ [El] , cosA[GJ]j [u] 


U! (D-67) 


It is now possible to define the elasticity matrix [S^] 
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[C] = 


1 

2 


E^i hh 


1 1 


2 E 2 I 2 


0 


[D] = 


1 1 
2 E (L-1) I (L-1) 


E L I L 


e l i l 


e l i l 


1 1 

2 Vl 


1 

R QUALITY 

1 1 
2 ¥r 

i i 

0 

Vr 2E (R+1) I (R+1) 


1 1 

1 

Vr e (r+i) i (r+i) 

’ ' ' " 2E . I 
1 1 

i l 

1 

Vr e (r+i) i (r+D 

2E I 
n n 


TABLE D-III 
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( 2 

[EI] - cosA [GJ]j [u] 

(D-68) 

Then 

!“si “ CS 2 3 1^1 ®- 69) 

The S2^ element of the [S^] matrix represents the angle of 
attack change in radians at station i due to the structural deflection 
of the wing caused by a unit loading at station J . 
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D.4 Wing Twist, 


The twist in a flexible wing can have several contributions or 
twists which can be divided into two classes: I) those which would be 
present even if the wing were rigid, and II) those due to inertia effects, 
thrusts or drag, and section pitching movements on the flexible wing. 


la I = a + 'a 


g 1 


St 


I gl 


(D-70) 


Class I: Aerodynamic Twists Ja 


8 ^ 


a) Built-in geometric twist due to camber or construction 
(including eventual dihedral contribution) or both. 

b) Interference twist (not considered in this study). 

c) Twist due to control surfaces deflection (flap, aileron, 
spoilers) . 

d) Apparent twist due to airplane rolling and/or pitching 
velocities (the angles of attack due to pitching velocities 
should be measured at 3C/4), 


Class II: Structural Twist due to Wing Deflections caused by Aerodynamic 

Loading which are independent of Wing Lift Distribution la I 

I g n I 

a) Vertical acceleration upon dry wing and internal fuel dead 
weights. 

b) Effects of airplane rolling and/or pitching acceleration 
upon dry wing dead weight, wing internal fuel dead weight. 

c) Section pitching moment with control surfaces in neutral 
position (C of the airfoil). 
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d) Incremental section pitching moment due to control surfaces 
deflection. 

The type of twist due to the effects of wing deflections arising 
from loads which are independent of wing angle of attack, such as those 
listed under Class II, may be computed with the aid of equation (D-63) 


a 


’ll 


2h 

cosA 


sinA [El] ] M | + cosA [GJ] 


(D-71) 


where M and T are the wing bending moments and torsion along the 
wing elastic axis due to the loadings of Class II. For cases a) and b) 
the knowledge of mass for each section of span 2h and the position of 
the corresponding C.M. w.r.t, the elastic axis must be known. 

No attempt is made here to compute the section pitching moment 
with control surfaces in neutral position. 


D.4.1 Twist due to Control Surface Deflection. 

The effect of a control surface deflection is of introducing an 
aerodynamic as well as a structural twist. 

The aerodynamic twist can be computed from the following expression 
for the lift produced by control deflection [Ref. 25]. 


where 


! C ii6 ~ 2 


(it -~e 0 -+-sine 0 ) 


0q = sin ^ 


yjU 1 - 5 ) 


(D-72) 
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0 O = cos -1 (2?- 1) 



c 


f 


control surface chord 



control surface deflection in radians 


Equation (D-71) applies for a 2-D lift curve slope of 2 tt . For a 
2-D lift curve slope of m^ , this equation can be written as 


L = ! 

V 


[rr-V sln6 ° )1 

16 

o v 


A 17 


Substituting the value for 0^ results in 


(D- 73) 


^ 16 


^ 1 - i~ccTs‘ 1 '(2l '■l) — h~~j_2 jyjui - d ) 


6 (D- 74) 


And dividing both sides by 


m. 


'a 


( i-i cos“ 1 (25- ii? ) 


(D-75 ) 


Structural twist derives from the section pitching movement due to 
control deflection. Once the spanwise bending moments and torsion 
distributions are known, equation (D-71) can be used to compute the 
section structural twist where M and T are replaced by the corre- 
sponding ones due to controls, M c and . 

In Figure D- 9 the symbols used are shown. 
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Figure D-9 - Structural Moments due to Control. 



C 


C 

m 0i 

M 

c, 

X 


T 

c. 

1 


Station pitching moment due to controls 


M _ = q 2h C C 

cO. m n . 

1 Oi 


Local chord length 

Local pitching moment coefficient 




sin A 
cosA 


The station pitching moments can be computed similarly to what was done 
in equations (D- 42) for the left wing, and (D-50) for the right wing. 

For stations 1 and n 


M 


cO 


1 


q h 


C 

m 01 


(D-76) 



n 


q h 



(D-77) 
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For stations 2 and n - 1 


M c 0 2 = q 2h 


M 


cO 


(n-1) 


= q 2h 


2 12 

c: c + -r tit c 

1 m Ql 2 1 m Q1 


12 2 
r C. n C + C C 

2 (n-1) m 0 ( n - ]_) n m 0n 


and so forth. 

And in general, for the left wing 



' i-l 

_ 

Kl = q 2h 

E 

. k=l 

2 12 
c c + - c. c 

k m Qk 2 1 m 0 i_ 


for the right wing 



' i+l 

K)* = q 2h 

V s 2 12 

2_J c, c + - c. c 

, k m ni 2 i m n , 

Lk=n Ok Oi 


Equations (D-80) and (D-81) are derived for the general case 
the control surface can be extended from tip to tip. 

In matrix form, (D-80) and (D-81) become 

= 2h q [Dll] 


) M I 

! cO I 


c 2 


"mO 


where 


[Dir> 


1/2 


1 

1 



1 

1/2 

1 



1 

1 

1/2 | 



* 4 

1 

4 4 4 

_ 1 _ 

1 .... 1/2 1 

1 

1 

1/2 

1 1 .... 1 



1 

1 

1 


1/2 1 .... 1 



I 

1 


1/2 1 



1 

t 


1/2 


(D-78) 

(D-79) 

(D-70) 

(D-81) 

where 

(D-82) 

(D-83) 
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If wind tunnel data are not available, theoretical expressions for 

C „ in terms of the deflection 6 may be used. 
mU 

From reference [25] 


“mO 1 6 ” " 2 


(sin0]'-~'Ysin20 ) 
U l " 


(D-84) 


Substituting sin0g and sin20g = 2sin0^ cos0q results in 


mO 1 6 


2 V 5(1- ?) - 2~V5(1--P<25- 1)' 


= - 2 


VYd- 


(D-85) 


Therefore, the structural bending moment and torsion are 


M c j = - 4 q h sinA [D I 1] 


Via'- S) 3 ' 




(D-86) 


1 T c I “ ~ 4 q h cosA [Dll] 


< 


5(1'- 5)' 


(D-87) 


and the structural twist is 


a 


s n’ 

. 2 , 


2h 

cosA 


sinA [El] [ M j + cosA [GJ] ]T 


I T ' 


8qh 2 [El] + cosA [GJ] 


[Dll] 


1(1"- 5); 




(D-88) 


The total twist due to control deflection is 



(D-89) 


a 


= a 


! +|a ! 


g 


II' 



cos ~ 1 ( 2 § 2=11 

TT 


•+ 


2 Vs(l- 5)‘ \ 


8 qh 2 


2 

sin A 
cosA 


[El] 





1 

+ cosA [GJ] 

[D I 1] 

tIui-JZ) 3 ' 


c/l 





1 ■- 



Equation (D-89) can be applied to both aileron and/or flap systems 
with no requirements for symmetry or antisymmetry in their displacement 
since, once [£}is specified, the deflection vector {6 ] can assume any 
configuration. 


D.4.2 Apparent Twist due to Rolling and Pitching Velocities. 

When the wing is at a skew angle, the apparent twist is computed 
similarly for both the rolling and pitching case, the only difference 
being in the way the station arms are computed. 

A plan view of the geometry of the quantity used in computing the 
apparent twist is given in Figure D-10. 

Consequently, the apparent twist is given by: 

1) for rate of'roll p 

(a ). = £ x. 

\ pji v 1 

2) for rate of pitch q 

K)i ■ - v (y 3 /4 )i 
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X^ = ” Y i tan A + f(y.) cos A 


63/4 A 


c + d cos A + X. 
P i 

C. 

— x 

y i " T 


x = d sinA + Y. 

x 


Figure D-10 - Planview of the Geometric Parameters Used in Computing 
the Apparent Twist due to Angular Velocities. 
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and in matrix form 





q l y 3/4| 


D.5 Unsymmetrical Flight Conditions. 

Because of the oblique wing, whenever an unsymmetric spanwise lift 
distribution occurs, the aircraft will experience accelerations about 
its rolling as well as pitching axis. 

A number of unsymmetrical flight conditions are usually investigated 
in structural design; for the purpose of this study, such investigation 
will be restricted to those which arise through the use of control 
surfaces on the wing, such as ailerons. 


The load distribution on an elastic wing associated with control 
deflections may be thought of as the summation of distributions from the 
following specific loadings. 


1) Symmetrical loading with controls in meutral position. 

2) Incremental loading due to controls deflection. 

3) Incremental loading associated with constant rolling and/or 
pitching velocity with controls in neutral position. 

4) Incremental loading caused by rolling and/or pitching angular 
acceleration. This loading results from the structural twist 
described under "Structural Twist" (class II b). 


In a steady roll and/or pitch condition, the span load distribution 
for the elastic wing is given by the first three loadings enumerated. 

For this condition the j Ct { values vary linearly and antisymmetrically 
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across the span from - (p + q tanA) at t ^ e ^- e ^ t (forward) tip, and 
b 

(p + q tanA) ^ at the other. 

The wing being at a skew angle, any unsymmetry in the spanwise lift 
distribution produces both pitching and rolling moments and, therefore, 
it introduces angular accelerations about the corresponding axes. Such 
a condition occurs when deflecting the desired wing control surface, 
unless some corrective action is taken in compensating the pitching 
moment by mean of the horizontal tail. 

We shall now investigate this case which implies non-zero p and/or 
q and divide it into three parts: initiation, steady state, and termi- 
nation of the motion. The second part, the steady roll and/or pitch 
condition, has already been discussed. The first and third ones differ 
only for the initial condition; no angular velocities for the initiation 
and a steady angular velocity for the termination. Therefore the analysis 
is the same for both cases. 

Let us therefore consider the initiation of the motion due to an 
instantaneous deflection of the ailerons. 

There will be contributions to the unsymmetric loading from all 
of the four loadings enumerated above. The first three ones having 
already been discussed, so only the loading due to the angular accelerations 
needs to be analyzed. The result of the angular accelerations is a new 
contribution to the structural twist due to inertia bending and torsional 
moments. The resulting twist distribution will superimpose an unsymmetric 
lift distribution on the already existing spanwise lift distribution. 

Because of the linear theory assumption this structural twist contribution 
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is independent from the othej* ones and, therefore, can be analyzed 
separately. 

Figure D-ll shows the geometry of the wing section center of mass. 

th 

The acceleration acting on the i center of mass is 

= x^ p + ynu q (D-90) 

and the corresponding loading 

F i = m i (x i p + ynu q) (D-91) 


Similarly to what done in "Structures Fundamentals" it is now 
possible to determine the pitching and rolling moments distributions and 
therefore the bending moment and torsion distributions. 


For the left wing 



M 

y 


1 


hfl 

2 2 


h m l 


2f(V + ^ *> 


m 


1 — — 

— (x 1 p + ym x q) (xrt^ + - tanA) 


(D-92) 


(D-93) 


and in general (for l<i£L where L = last station to the left of air- 


craft centerline) 
i-1 


M 


= 2h m^k - i)(x k p + ym k q) - m jL (x i p + ynu q) (D-94) 


i k=l 


M 


i-1 

= £-*>. 


i K=1 


i [^ki ” (e k " Xm k> + e i ](x k 


(D-95) 


m. 


- y ^i + 2 tan ^( x i P + yn^ q) 
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(C.M.) ± 

xni. 

i 

*“i 

X. 


Center of Mass of i wing section 

distance of (C.M.). from elastic axis 
1 

c + d cosA + X. - (e, -xm.) 
p x v 1 1 

- Y. tanA + f (y. ) 
x y x 


Figure D-ll - Wing Section Center of Mass, 
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For the right wing 


M * * 2 I ( *n P + >™ n 
n 


(D-96) 


M ■ = 
y n 


Y ( xm n - 2 tan A) (x n p + yn^ q) 


(D-97) 


and in general for (n>iSR) where R - first station to the right of 

aircraft centerline 
i+1 


M x = X) 2h m^k - 1)0^ P + q) + jr p + ynu q) (D-98) 

i k=n ,v ’ ' r x 


i+1 

M y = 2 - • e i) + ^ k P + ^ 

J x k=n 


J=I1CN 

1 

0*1" 

1 

tanA) (x i 

P + 

ym i q) 

and in matrix form 





l M 1 = 
lx! 

I* 1 ! 

P + 

1 A 2 

1 * 

l M r 

1 y) 

[b ij 

p + 

}b 2 

u 


where 


(a l' 
( A I, 

- h[A][-m f J 

t* 

1 

1 


1 A 2 ^ 

r l i 

= h[A]['m £ J 

W 


|B1| 


*1 



i B2 l 

= [B]rm f J 

jym 

l 

1 



= mass matrix (diagonal). 


[B] < 


[E ] 

; co] 

i 



[0] 

i [f] . 


(D-99) 


(D-100) 


(D-101) 


where [E] and [F] are reported in Table D- IV . 
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TABLE D - IV 


ORIGINAL PAGE IS 
OP POOR QUALITY 



Since 


iv 


= ! M J cosA + ! M ! sinA 


x : 


^ T J = !M ! cos A - }M ! sinA 
> ! y ) ! x ! 

it is now possible to compute ^ Q' ] and jo; | from equation (D-71). 

.. i s n )^ 

The solution of the following equation will give the lift distribu- 
tion arising from the twist contribution due to p and q 


[S 1 ] - [S 2 ] 



+ 



g II q 


(D-103) 


It is now possible to compute the aerodynamic pitching and rolling 
moment due to p and q by integrating over the entire span the lift 
force times the corresponding moment arm. 


M x = h |_1 2 2 ... 2 lj r-x-J \ &\ 
' j 

M = h [I 2 2 ... lj ["y-J \& | 

j J 

where j = p , q . 


(D-104) 


By assuming now that p and q are small quantities, it is 
possible to linearize even the p and q contributions to aerodynamic 
pitching and rolling moments 


M 

y 


3M 


L = 


x 


3p 


P + 



3M 



P + 



(D-105) 
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where 


M 

3M x. 



and 


j = P , 4 
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